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EXECUTIVE  SUMMARY 


Compared  to  many  other  forging  methods,  radial  forging  is  relatively 
new  having  been  commercialized  late  in  the  twentieth  century.  Radial 
forging  may  be  characterized  as  a  multi-hammer,  open  die  forging  process  for 
producing  solid  and  hollow,  rotational  and  non-rotational,  axisymmetric 
parts.  Two  types  of  machines  are  in  common  usage  employing  either  two  or 
four  hammers.  Of  the  various  radial  forging  machines,  Gesellschaft  fuer 
Fertigungstechnik  und  Maschinenbau  GmbH  (GFM)  is  the  most  widely 
recognized  manufacturer. 


This  report  summarizes  the  state-of-technology  in  radial  forging.  A 
brief  presentation  of  the  computer  program,  RFORGE,  which  was  developed 
by  Battelle  under  an  earlier  contract  is  included  (RFORGE  User’s  Manual  is 
included  in  Appendix  A)  as  are  initial  results  of  a  reduced  scale,  isothermal 
simulation  of  the  radial  forging  process  using  the  present  preform  design. 
The  simulation  was  run  using  the  computer  program  DEFORM  which  was 
also  developed  by  Battelle  specifically  for  simulating  metal  forming  processes. 
RFORGE  was  found  to  be  useful  for  situations  where  detailed  process  analyses 
are  not  required.  However,  RFORGE’s  lack  of  post-processing  capabilities 
requires  that  the  user  manually  process  and  interpret  the  results.  On  the  other  ^ 
hand,  DEFORM  was  found  to  be  very  useful  in  modeling  the  radial  forging 
process  and  offers  the  user  detailed  post-processing  and  information  on 
strain,  strain  rate,  temperature,  and  metal  flow.  To  simplify  the  finite  element 
process  model  used  in  this  study,  the  simulation  was  approximated  as  an 
isothermal,  axisymmetric  process.  Simulation  results  are  included. 


Data  supplied  by  Benet  Laboratory  was  analyzed  using  statistical - - 

methods  to  determine  what  significant  relations  exist  between  the  forged  and  q* 
heat  treated  gun  tube  and  resultant  mechanical  properties.  Due  to  probable  □ 

confounding  effects  of  the  forging  and  heat  treatment  processes  no  significant _ 

relations  were  found.  It  is  recommended  that  a  more  detailed  experimental  _ _ 

study  be  carried  out  which  isolates  the  forging  process  from  post-forging  heat  Jfc — . 
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treatment  in  order  to  gain  a  better  understanding  of  how  the  radial  forging 
process  affects  the  mechanical  properties  of  the  forged  gun  tube.  A  simple  true 
strain  method  of  calculating  forging  reduction  was  used  to  analyze  the  plastic 
strains  at  the  breech  end  to  develop  a  new  preform  design  which  is  expected 
to  yield  uniform  amounts  of  strain  at  each  point  in  the  gun  tube  after  radial 
forging. 
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INTRODUCTION 


Metal  forming  processes  may  be  divided  into  two  broad  classifications 
depending  on  their  flexibility.  Lange  has  proposed  such  a  scheme  where 
metal  forming  is  divided  into  free-flow  forming  and  restricted-flow  forming 
processes  (Reference  1).  According  to  this  definition,  the  final  shape  and 
tolerance  of  the  workpieces  in  free-flow  forming  is  dependent  only  on  the 
relative  motions  and  positions  of  the  dies.  Restricted-flow  forming  on  the 
other  hand  is  independent  of  machine  motion  and  is  a  function  of  the  die 
impression  only  based  on  Lange's  classification.  Based  on  these  definitions  it 
is  readily  seen  that  radial  forging  may  be  broadly  classified  as  a  free-flow 
forming  process.  Additionally,  it  is  apparent  that  radial  forging  is  a  flexible 
forming  process  in  that  a  variety  of  different  part  geometries  can  be  produced 
without  any  requirements  for  die  change. 

Radial  forging  often  referred  to  as,  and  confused  with,  rotary  forging  is 
a  relatively  new  forging  process.  The  term  rotary  forging  is  used  due  to  the 
fact  that  in  some  operations  the  workpiece  is  rotated  during  the  forging 
operation.  The  term,  radial  forging,  appears  to  be  widely  used  in  Europe  to 
describe  the  process  while  both  radial  and  rotary  forging  are  commonly  used 
in  the  United  States  to  describe  the  process.  The  process  is  depicted  in  Figure  1 
below. 

A  similar  process  that  may  be  classified  as  a  form  of  radial  forging  is 
rotary  swaging.  Swaging  differs  from  radial  forging  in  that  the  dies  are 
mounted  in  a  headstock  which  holds  the  dies  and  a  set  of  planetary  rollers. 
Typically,  two  or  four  dies  are  used.  The  end  of  the  dies  which  contacts  the 
rollers  is  shaped  as  a  cam.  As  the  headstock  rotates  at  a  high  frequency  the 
dies  are  activated  by  the  rollers  and  returned  by  centrifugal  force.  The  die 
travel  is  changed  by  adding  or  removing  shim  stock  between  the  dies  and 
rollers. 


1 


Radial  forging  in  its  present  form  is  essentially  an  open  die  forging 
process  using  sets  of  counterblow  hammers  or  dies  where  tubular  (solid  or 
hollow)  or  square  billets  are  reduced  in  cross-section  to  straight  or  tapered 
forms.  For  processing  of  tubular  billets,  a  mandrel  is  required.  Radial  forging 
is  commonly  done  using  2,  3,  or  4  counteracting  dies  located  in  one  plane. 
The  simplest  and  most  common  orientation  is  the  vertical  plane.  Most 
forging  applications  for  radial  forging  involve  simple  reduction  of  the  cross- 
section  of  axisymmetric  parts.  However,  increased  use  of  the  process  is  being 
made  in  small  lot  production  of  complex  "T",  "X",  and  "Y”  shapes.  Radial 
forging  is  also  capable  in  some  cases  of  producing  parts  of  non-symmetric 
cross-section.  However,  the  machine  must  be  modified  to  withstand  the 
torque  and  moment  arms  developed  by  the  out-of-plane  forging  conditions 
(Ref.  2).  Additionally,  this  necessitates  that  the  machine  be  capable  of 
independent  adjustment  of  the  individual  hammers. 

Rotary  forging  was  developed  from  open  die  forging  where  solid  billets 
are  reduced  between  a  flat  upper  die  and  a  "U"  or  "V"  shaped  lower  die  in  a 
mechanical  or  hydraulic  press.  An  overhead  crane  or  track  mounted 
manipulator  is  then  used  to  rotate  and  feed  the  billet  between  subsequent 
press  hits.  This  process  is  depicted  below  in  Figure  2: 
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Figure  2.  Open  Die  Forging  Process  (Ref.  3). 


However,  there  are  a  number  of  problems  inherent  in  the  open  die 
process  depicted  in  Figure  2.  The  process  is  quite  slow  due  to  the  limitation  of 
stroking  rate  of  the  press  being  used.  Due  to  the  die  configuration,  edge  and 
surface  cracks  or  bursts  are  likely  to  develop  due  to  the  tensile  stresses  that  can 
develop  in  the  unsupported  edges  of  the  workpiece.  Additionally,  significant 
heat  loss  from  the  billet  to  the  dies  may  require  that  multiple  reheats  be 
performed  to  prevent  the  material  flow  stress  from  becoming  too  large. 
Center  bursting  and  surface  cracking  problems  may  be  minimized  by  using 
"V"  shaped  dies  but  this  limits  the  amount  of  reduction  available  in  one  pass 
and  may  require  frequent  die  changes  between  passes  (Ref.  4).  Many  of  these 
problems  have  been  overcome  by  the  present  radial  forging  machines. 

The  basic  configuration  of  the  radial  forging  machine  was  initially 
developed  by  Dr.  Bruno  Kralowetz  who  founded  GFM  in  Steyr,  Austria  in 
1946  (Ref.  5).  The  first  four  hammer  working  model  of  the  machine  was 
introduced  in  1960.  GFM  has  since  installed  over  600  machines  worldwide.  A 
cross-sectional  view  of  the  forging  box  of  a  GFM  machine  is  shown  below  in 
Figure  3. 
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Figure  3.  Cross-section  of  the  GFM  Radial  Forging  Machine  (Ref.  6). 


As  can  be  seen,  the  basic  premise  of  the  machine  is  that  four  dies 
attached  to  hammers  simultaneously  press  the  workpiece.  It  is  important  to 
note  the  distinction  that  the  dies  actually  press  the  workpiece  although  it 
would  appear  that  with  the  high  stroking  rate  that  the  dies  actually  act  like 
hammers.  However,  the  deformation  rates  encountered  in  radial  forging  are 
within  those  typically  encountered  for  press  operations.  The  radial  forging 
process  has  a  high  stroke  rate  due  to  the  fact  that  the  hammers  are  spaced  very 
close  to  the  workpiece,  requiring  only  a  short  period  of  time  between 
subsequent  strokes.  As  the  radial  forging  machine  is  simply  a  short-stroke 
mechanical  press,  the  die  velocity  at  any  point  in  the  stroke  may  be  calculated 
by  using  equation  (1)  below  (Ref.  7).  Additionally,  many  machines  have  an 
automated  stock  handling  system  incorporated  into  the  machine  to  provide 
coordinated  rotation  and  feed  of  the  stock  between  hits.  Thus,  only  a  small 
portion  of  the  workpiece  undergoes  deformation  at  a  given  time.  The 
advantage  of  using  four  hammers  is  that  the  probability  of  edge  cracks 
developing  from  tensile  stresses  is  quite  small  as  the  surface  is  under 
compression.  Due  to  the  counteracting  motion  of  the  hammer  dies,  there  are 
virtually  no  stresses  being  transmitted  to  the  forging  machine  frame  and 
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since  only  a  small  portion  of  the  workpiece  is  being  deformed  at  one  time,  the 
forging  forces  required  will  be  relatively  low.  Heat  loss  from  the  workpiece  is 
minimized  and  in  many  cases,  with  proper  selection  of  process  parameters, 
the  temperature  may  increase  from  deformation  as  the  hammers  contact  the 
workpiece  typically  from  150-200  times  per  minute.  This  is  advantageous  as  it 
eliminates  the  need  to  reheat  the  billet  between  passes.  The  need  for  die 
changes  are  also  eliminated  with  the  radial  forging  machine  as  the  amount  of 
reduction  may  be  varied  merely  by  changing  the  stroke  length  of  the 
hammers.  On  many  machines  this  may  be  done  by  simple  reprogramming  of 
the  machine.  Another  advantage  of  the  radial  forging  process  is  that  the 
internal  soundness  of  cast  billets  is  improved  and  internal  voids  closed. 


(1) 


V~ 


nNW  S__{ 
30  VW 


(Ref.  10) 


where:  W  =  Ram  Location  from  Bottom  Dead  Center  (inches) 

S  =  Stroke  Length  (inches) 

N  *  Number  of  Strokes  per  Minute 


Present  day  rotary  forging  equipment  is  best  characterized  by  the 
number  of  hammers  used  in  the  machine.  Two  and  four  hammer  machines 
being  most  common.  Typical  hammer  arrangements  for  each  machine  are 
shown  below  in  Figures  4,  5,  and  6.  Another  classification  scheme  would  be  to 
group  machines  on  the  basis  of  whether  the  hammers  are  arranged  in  the 
vertical  or  horizontal  plane.  However,  radial  forging  machines  of  the  latter 
type  are  not  common. 


Variables  in  Radial  Forging 

Shown  below  in  Table  1  is  a  listing  of  the  major  variables  found  in 
radial  forging  which  must  be  controlled  to  obtain  desired  finish  product 
properties. 
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Table  1.  Significant  Variables  in  Radial  Forging  Process. 

1.  Forging  Stock  or  Pre-form  Variables: 

-  Material  Flow  Stress 

-  Microstructure 

-  Stock  Dimensions 

-  Material  Physical  Properties  (e.g.  thermal  properties) 

-  Scale  on  Surface  of  Stock 

-  Temperature  of  Workpiece 

2.  Machine  and  Tooling  Variables 

-  Die  Geometry 

-  Die  Material  and  Properties 

-  Mandrel  Properties  and  Design 

-  Use  of  Lubricant  on  Mandrel 

-  Stroking  Rate 

-  Length  of  Stroke 

-  Machine  Tonnage  Capacity 

3.  Process  Variables 

-  Billet  Temperature 

-  Die  Stroke  Velocity/Contact  Time  on  Workpiece 

-  Percent  Reduction  of  Workpiece 

-  Temperatures  Generated  in  Work  Zone 

-  Feed  Rates  (Axial  and  Longitudinal) 
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RADIAL  FORGING  MACHINES 


Two  Die  Radial  Forging  Machines 

Two  die  radial  forging  machines  are  quite  similar  to  traditional  open 
die  forging  presses  in  that  two  dies  are  used.  Radial  forging  machines  on  the 
other  hand  have  a  much  higher  stroking  rate  than  traditional  mechanical 
presses. 

I.  Fenn  Manufacturing  Radial  Forging  Machine 

A  schematic  of  a  two  hammer  radial  forging  machine  built  by  Fenn 
Manufacturing  Company  under  license  from  Usine  de  Wecker  is  shown 
below  in  Figure  4  (Ref.  8).  As  of  12/31/89,  the  manufacture  of  these  machines 
in  Europe  was  transferred  to  GFU-Maschinenbau  GmBH  of  West  Germany 
(Ref.  8).  With  suitable  alteration,  the  machine  can  also  operate  with  four  dies. 
As  can  be  seen,  the  dies  are  mounted  on  rocker  arms  driven  by  a  single 
driveshaft  via  eccentric  cams.  The  dies  are  returned  by  spring  action  after 
deforming  the  workpiece.  The  distance  between  the  dies  is  controlled  by 
adjusting  a  worm  type  bushing  which  controls  the  throw  of  the  dies  or 
chuckheads  which  also  advance  the  part  in  the  axial  direction.  Due  to  the 
rapid  stroking  of  the  dies,  significant  heat  buildup  can  occur  in  the  dies.  To 
avoid  overheating,  the  dies  are  water-cooled.  As  can  be  seen  in  Figure  4,  the 
Fenn  Machine  is  a  relatively  simple  mechanism.  Maximum  operating  speed 
is  1000  strokes  per  minute  with  forging  loads  of  15,000  KN  per  tool  (Ref.  8). 

II.  Ruthner  Radial  Forging  Machine 

The  two  die  hydraulic  radial  forging  machine  from  the  Andritz- 
Ruthner  Industrieanlagen-Akliengeselleschaft  of  Austria  was  first  introduced 
in  1981  (Ref.  9).  The  forging  box  consists  of  two  horizontally  arranged 
hydraulic  cylinders  upon  which  the  forging  dies  are  mounted.  Some  of  the 
machine  features  include  1,250  tons  capability  per  tool,  automatic  tool¬ 
changing  and  pre-heating,  and  a  hydrostatic  mounting  system  for  the 
manipulators.  The  uniqueness  of  the  hydrostatic  mounting  system  stems 


from  the  fact  that  the  base  of  the  manipulators  ride  entirely  upon  a  layer  of  oil 
in  order  to  eliminate  all  metal  to  metal  contact  and  wear  (Ref.  9). 


Figure  4.  Two  Die  Radial  Forging  Machine  (Ref.  8). 


Figure  5.  Three  Die  Radial  Forging  Machine  (Ref.  4). 
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Due  to  the  two  hammer  layout  of  the  Ruthner  machine,  the  forging  of 
round  stock  requires  multiple  forging  passes.  Initially  the  billet  is  pass  forged 
to  achieve  an  octagon  shape  slightly  larger  than  the  finish  forged  dimension 
while  the  final  pass  forges  the  finish  dimension  using  short  strokes  to  obtain 
small  depth  penetration  at  the  surface  (Ref.  9). 

Four  Die  Radial  Forging  Machines 

In  general,  the  four  die  radial  forging  machines  are  more  My  used 
than  two  die  machines.  One  of  the  major  advantages  of  the  four  achines 
is  that  square  and  round  shapes  may  be  produced  in  one  pass  wheit^  the  two 
die  machines  often  require  multiple  forging  passes.  In  general,  four  die 
machines  are  more  complex  and  costly  than  two  die  machines. 

I.  Kocks  Four  Die  Swing  Forging  Machine 

This  machine  differs  from  other  radial  forging  machines  in  that  the 
dies  are  activated  by  mechanical  linkages  and  eccentric  cams  as  shown  in 
Figure  5.  Most  four  die  machines  are  driven  by  either  hydraulic  pistons  or 
mechanical  connecting  rods.  Additionally,  the  Kocks  machine  differs  from 
other  radial  forging  machines  in  that  it  is  intended  mainly  for  hot  reduction 
of  square  or  rectangular  billets  on  a  continuous  basis,  similar  to  rolling  (Ref. 
4). 


Each  die  is  operated  independently  by  three  linkages  mounted  on 
eccentric  cams.  The  dies  are  operated  in  pairs  with  the  second  pair  operated  at 
180  degrees  to  the  first  pair.  The  machine  is  self- feeding  and  requires  no 
manipulators  to  feed  the  stock  as  the  dies  pull  the  material  through  the  forge 
zone  after  deformation.  The  dies  initially  sink  the  material  as  shown  in 
Figure  6a.  With  continued  travel  of  the  dies,  the  material  is  finish  forged  and 
pulled  through  at  the  end  of  the  die  rotation. 
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Figure  6.  Kocks  Swing  Forging  Machine  (Ref.  4). 

II.  GFM  Radial  Forging  Machine 

The  four  die  configuration  of  the  radial  forging  machine  is  probably  the 
most  common  machine  in  use.  Of  these,  GFM  is  best  known  for  radial  forging 
machines. 

The  first  radial  forging  machines  built  by  GFM  were  based  on  a 
pneumatic-hydraulic  drive  system  which  closely  approximated  the  action  of  a 
counterblow  hammer.  A  schematic  of  the  machine  is  shown  below  in  Figure 
7.  The  dies  were  activated  by  a  high  pressure  nitrogen  system  connected  to 
individual  pistons.  This  enabled  the  dies  to  achieve  high  impact  energies  and 
velocities  up  to  30  feet/second  (Ref.  10).  Synchronization  was  provided  by  a 
hydraulic  system  to  maintain  die  travel.  Upon  completion  of  the  forming 
stroke,  the  dies  were  retracted  by  the  hydraulic  system  which  also  acted  to 
compress  the  nitrogen  gas  (Ref.  10).  The  machine  was  capable  of  preheating 
the  stock  with  the  addition  of  an  induction  heater  mounted  at  the  entrance  to 
the  forge  box.  This  system  has  been  replaced  in  favor  of  the  mechanical  drive 
system  described  below. 

The  present  four  hammer  GFM  radial  forging  machine  is  essentially  a 
mechanically  driven  short  stroke  press.  Due  to  the  close  proximity  of  the 
hammers  in  relation  to  one  another,  very  high  stroking  rates  are  achieved. 
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This  is  a  major  change  from  the  GFM  machine  described  above.  Each  die  is 
mounted  on  a  connecting  rod  which  is  driven  by  an  adjustable  eccentric  shaft 
which  enables  variable  stroke  lengths.  The  eccentric  shafts  are  driven  via 
gears  which  synchronize  the  amount  of  die  travel.  Both  round  and 
rectangular  cross-sections  may  be  reduced  as  die  travel  may  be  adjusted  either 
in  sets  or  in  unison  via  gears  in  the  forging  box.  GFM’s  largest  machine  used 
in  the  West,  the  SX-65,  is  capable  of  125  strokes  per  minute  with  a  forging  load 
of  25,000  KN  per  tool  (Ref.  6).  A  larger  model,  the  SX-85,  has  been  delivered  to 
the  Soviet  Union. 


The  construction  of  the  GFM  radial  forging  machine  is  radically 
different  from  that  of  typical  open-die  forging  presses.  The  forging  action  of 
the  radial  forge  takes  place  within  a  vertically  arranged  forging  box  which 
houses  the  four  hammers  and  drives  at  right  angles.  Due  to  the  opposing 
motion  of  the  hammers,  no  forces  are  transmitted  to  the  machine  itself. 
Therefore  the  machine  foundation  remains  virtually  free  from  stress  and 
vibration  (Ref.  6). 


Figure  7.  GFM  Four  Die  Radial  Forging  Machine  (Ref.  10). 


11 


III.  SMS  Hasenclever  Radial  Forging  Machine 

A  rather  unique  four  hammer  radial  forging  machine  was  developed 
by  SMS  Hasenclever  in  collaboration  with  the  Technical  University  of 
Stuttgart  to  forge  tubular  and  solid  billets.  A  cross-sectional  view  of  the 
RUMX  2C00  is  shown  below  in  Figure  8.  The  uniqueness  of  the  RUMX  2000 
derives  from  the  fact  that  it  allows  individual  control  of  each  ram  stroke  and 
adjustment  of  the  stroke  length  and  position  between  each  stroke.  Another 
feature  of  the  RUMX  2000  is  a  programmable  tool  change  system  which 
enables  quick  change  capability  of  the  dies  during  the  process  which  allows  a 
workpiece  of  550  square  mm  to  be  reduced  to  65  mm  diameter  without 
interruption  (Ref.  1).  With  appropriate  modifications,  the  RUMX  2000  is  also 
capable  of  forging  non-symmetric  workpieces. 

The  hammers  are  hydraulically  driven  and  have  a  maximum  capacity 
of  10  MN.  Machine  control  is  provided  by  a  DEC  PDP  11/73  computer  (Ref.  1). 
The  control  computer  also  generates  process  plans  and  appropriate  processing 
control  program  from  the  finished  workpiece  geometry  using  the  computer 
program  "PRORUM"  as  described  below. 


Radial  Forging  Machine  Drives 

A  review  of  the  common  radial  forging  machines  show  that  two  drives 
are  in  common  use:  mechanical  and  hydraulic.  Of  these  two,  mechanical 
drives  are  predominant.  The  primary  reason  for  the  popularity  of  the 
mechanical  drive  is  due  to  the  inherently  faster  speed  of  the  drive  compared 
to  a  hydraulic  system.  However  hydraulic  drives  offer  greater  flexibility  in 
forging  prismatic  and  non-symmetric  cross-sections  due  to  the  possibility  for 
individual  control  of  each  hammer's  stroke. 
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Figure  8.  SMS  Hasenclever  Four  Die  Radial  Forging  Machine  with 
Tool  and  Workpiece  Motions  (Ref.  11). 
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Radial  Forging  Automation 


One  common  feature  of  most  radial  forging  machines  is  integral, 
automatic  stock  handling  by  one  or  two  manipulators.  For  the  case  where  the 
workpiece  must  be  forged  over  its  entire  length  in  one  heat,  two  chuckheads 
are  normally  used  (Ref.  12).  In  the  case  of  two  manipulators,  one  manipulator 
serves  to  feed  the  workpiece  while  the  second  manipulator  merely  holds  the 
workpiece  and  maintains  centering.  The  stock  manipulators  serve  to  pro\§de 
both  longitudinal  and  rotational  movement  to  the  workpiece.  Workpiece 
movement  occurs  between  hammer  blows  to  prevent  twisting  of  the 
workpiece.  Unlike  open  die  forging,  the  manipulators  constantly  maintain 
the  workpiece  center  position  between  the  dies  irregardless  of  the  cross- 
section  reduction.  The  manipulators  perform  both  axial  and  rotational  feed 
between  strokes. 

Many  radial  forging  installations  have  completely  automated  processes 
integrated  into  flexible  manufacturing  systems  (FMS).  A  typical  example  of  an 
automated  radial  forging  operation  is  the  line  at  Carpenter  Technology’s  forge 
plant  in  Reading,  PA  which  was  installed  in  1983.  This  line  is  depicted  in 
Figure  9  below.  Another  example  of  an  automated  radial  forging  line  that 
integrates  forging,  heat  treatment,  and  automatic  stock  handling  operations, 
is  located  at  the  U.S.  Army's  Watervliet  Arsenal  in  Watervliet,  NY. 
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Figure  9.  A  Radial  Forging  Line  at  Carpenter  Technology's  Plant  (Ref.  13). 

Forging  Over  A  Mandrel 

The  primary  difference  between  forging  solid  and  hollow  workpieces  is 
that  a  mandrel  must  be  used.  For  hot  and  cold  forging  of  seamless  tubes  in 
radial  forging,  a  retractable  mandrel  must  be  employed  to  support  the  inner 
wall  of  the  workpiece.  An  additional  difference  is  that  three  deformation 
zones  exist  during  tube  forging  whereas  solid  forging  has  only  two 
deformation  zones.  This  is  depicted  below  in  Figures  10  and  11.  Typically  the 
mandrel  is  a  hardened  piece  of  tool  steel  with  internal  water  cooling  to 
prevent  seizure  due  to  exposure  to  the  high  temperatures  and  pressures 
generated  during  the  forging  process.  The  outer  diameter  of  the  mandrel 
corresponds  to  the  finished  inner  diameter  of  the  forged  workpiece.  Two 
types  of  mandrels  are  commonly  used  in  radial  forging.  The  first  is  a  "short" 
or  stationary  mandrel  which  is  used  to  forge  tubes  having  a  constant  internal 
bore  (Ref.  12).  The  second  type  is  termed  a  "long"  or  moveable  mandrel  and  is 
used  to  forge  tubes  having  stepped  inner  and  outer  diameters  (Ref.  12).  The 
long  mandrel  moves  in  progression  with  the  workpiece  through  the 
deformation  zone  while  the  short  mandrel  remains  stationary.  Due  to  the 
fine  surface  finishes  from  cold  forging,  some  work  has  been  done  in  finish 
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forging  the  rifling  in  small  caliber,  thin  wall  gun  barrels  (Ref.  14).  Due  to  the 
presence  of  high  temperatures  and  pressures  at  the  mandrel/ workpiece 
interface,  in  many  cases  adequate  lubrication  or  an  insulating  material  may 
need  to  be  applied  to  the  mandrel  to  prevent  sticking. 

The  stresses  on  the  mandrel  were  analyzed  by  Ragupathi  et  al.  (Ref.  16) 
using  the  computer  program,  ADINA,  which  is  based  on  the  finite  element 
method.  An  examination  of  the  stresses  on  the  mandrel  show  that  it 
subjected  to  both  thermal  and  mechanical  stresses.  Thermal  stresses  are  set-up 
due  to  the  creation  of  large  thermal  gradients  between  the  water-cooled 
mandrel  and  the  heated  workpiece.  Mechanical  stresses  arise  due  to  the 
intermittent  loading  from  the  hammers.  Ragupathi  et  al.  (Ref.  16)  reported 
that  the  thermal  stresses  remained  relatively  constant  during  the  process  due 
to  the  constant  nature  of  the  temperature  gradients  whereas  the  mechanical 
stresses  could  be  modeled  using  a  sinusoidal  waveform.  Since  the  dies 
operate  at  200  strokes /minute,  both  the  dies  and  mandrel  are  subject  to 
fatigue  failure  (Ref.  17). 
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Figure  10.  Deformation  Zones  in  Radial  Forging  of  Tubes  (Ref.  15). 
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Figure  11.  Deformation  Zones  in  Radial  Forging  of  Solids  (Ref.  15). 
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RADIAL  FORGING  RESEARCH 


The  first  significant  paper  was  written  by  Sachs  and  Baldwin  in  1946 
where  they  analyzed  the  stress  states  in  tube  drawing  (Ref.  18).  While  tube 
sinking  is  not  a  radial  forging  process  per  se,  it  is  quite  relevant  in  that  the 
material  deformation  and  stress  states  are  similar  to  those  in  tube  sinking  in 
radial  forging.  Sachs  and  Baldwin  were  the  first  to  apply  the  slab  approach  to 
analyzing  the  tube  sinking  process.  The  resulting  expression  for  the  draw 
stress  is  shown  below  in  equation  (2). 
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Their  result  shows  that  the  draw  stress  is  a  function  of  the  tube 
reduction  and  the  existing  friction  conditions  at  the  workpiece  and  die 
interface. 

The  next  noteworthy  paper  was  presented  by  Kegg  in  1964  on  the 
mechanics  of  rotary  tube  swaging  (Ref.  19).  Kegg  attempted  to  develop  a 
model  for  the  kinematics  of  the  heaastock  mechanism  and  for  predicting  the 
power  requirements  for  the  process  and  estimating  the  maximum  reduction. 
The  model  that  Kegg  used  in  his  analysis  for  determining  the  die  pressure 
assumed  that  the  radial  forging  process  was  similar  to  plane  strain  extrusion 
with  a  small  inlet  angle.  The  resulting  equation  from  the  slab  analysis  shows 
that  the  die  pressure  is  a  function  of  the  back  push  stress  material  flow  stress 
and  the  deformation,  x. 

(3)  <r,  =^<Tt+-^r<70j 

Examining  equation  (3)  above  it  can  be  seen  that  the  result  is  very  similar  to 
the  case  for  extrusion.  However,  this  analysis  ignores  the  effect  of  the 
finishing  zone  which  has  been  shown  to  be  a  significant  effect  in  radial 
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forging  (Ref.  15).  Additionally,  the  assumption  that  the  die  inlet  angle,  a,  is 
small  is  not  always  valid  in  radial  forging  practice. 

No  further  work  was  done  in  this  area  until  1974  when  researchers  at 
Battelle  Columbus  Laboratories  analyzed  the  radial  forging  process  for 
manufacturing  gun  barrels  for  the  United  States  Army  (Ref.  15).  This  effort 
was  the  first  to  extensively  mathematically  model  the  radial  forging  process. 
As  part  of  the  work,  slab  and  upper  bound  analysis  of  the  forging  of  gun  tubes 
was  performed  to  predict  the  stresses,  strains,  and  temperatures  during  the 
forging  process.  The  resultant  work  was  incorporated  into  a  FORTRAN  based 
computer  program  called  RFORGE  which  enabled  calculation  of  the  neutral 
plane,  forging  loads,  strain,  strain  rates,  and  temperatures  generated  in  the 
workpiece  during  the  forging  operation.  RFORGE  is  limited,  however,  to  the 
forging  of  tubular  workpieces. 

The  analysis  carried  out  by  Altan  et  al.  (Ref.  15)  assumed  that  the  radial 
forging  process  was  divided  into  three  distinct  zones:  sinking,  forging,  and 
sizing.  Additionally,  complete  contact  around  the  periphery  was  assumed 
between  the  dies  and  workpiece  such  that  no  flash  was  generated.  An 
important  departure  from  previous  analyses  was  that  the  coefficient  of 
friction  was  replaced  by  the  friction  factor.  Use  of  the  friction  factor  enables 
determination  of  the  friction  conditions  existing  in  metal  forming  by  the  use 
of  the  relatively  simple  ring  test.  Additionally,  the  ring  test  was  used  to 
generate  the  flow  curves  for  the  work  material  at  warm  and  hot  working 
temperatures. 

The  resultant  equations  obtained  for  the  sinking  zone  from  the  slab 
analysis  were  found  to  be  for  the  axial  stress  and  die  pressures  respectively 
(Ref.  15): 
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(5)  P^'Ma-a^icosa)2 
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The  resultant  equations  for  axial  and  die  pressures  in  the  forging  zone 
were  found  to  be  respectively  (Ref.  15): 

(AlMJLM] 

(Rt-Rm)(R  +  Rm) 

(7)  P  =  a-at 

The  resultant  equations  for  the  axial  and  die  pressures  in  the  sizing 
zone  were  found  respectively  to  be  (Ref.  15): 

(9) 

Lange  developed  a  computer  program  for  optimizing  the  radial  forging 
process  for  the  four  hammer  hydraulic  machine,  RUMX  2000,  designed  at  the 
University  of  Stuttgart  and  SMS  Hasenclever  (Ref.  1).  The  program  was  called 
"PRORUM"  and  using  the  geometry  of  the  final  machined  part,  strains, 
machining  allowance,  available  raw  material  dimensions,  and  available  tool 
sets,  the  inputs  for  an  automated  process  planning  were  generated.  PRORUM 
was  a  significant  breakthrough  in  the  application  of  NC  to  metal  forming  in 
that  all  necessary  forming  steps  in  the  part  program  were  generated  from  the 
input  data.  From  the  above  data,  PRORUM  adapted  the  geometry  to  the  radial 
forging  process,  determined  a  suitable  preform  size,  and  optimum 
deformation  sequence.  In  the  next  phase  of  the  program,  the  individual 
process  steps  are  determined  including  axial  and  angular  feeds,  ram  motions, 
and  which  of  the  available  tool  sets  should  be  used  to  obtain  the  desired 
degree  of  plastification  of  the  cross-section.  The  program  is  of  an  interactive 
nature  and  is  capable  of  displaying  the  workpiece  at  various  intermediate 
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steps.  Additionally  PRORUM  determines  if  there  is  any  danger  of  tool 
collision  based  on  the  tool  sets  chosen  for  the  process. 

Lange  further  extended  the  program  PRORUM  to  cover  optimization 
of  stepped  shaft  production  for  radial  forging  and  turning  combined  into  a 
complex  FMS  cell  (Ref.  20).  The  program,  "VORUM"  was  similar  to 
PRORUM  in  that  it  was  specifically  written  for  the  SMS  Hasenclever  RUMX 
2000  Forging  machine.  VORUM,  a  multi-function  program,  is  able  to  select 
the  processing  steps  for  a  complex  FMS  cell  which  allows  optimum  use  of 
both  the  radial  forging  machine  and  the  lathe  to  produce  a  given  stepped 
shaft.  Based  on  the  overhead  and  production  costs  input  as  data  to  the 
program,  VORUM  determines  all  possible  process  combinations  based  on  a 
breakdown  of  the  part  geometry  and  then  selects  the  sequence  of  forging 
and/or  turning  that  offers  minimum  processing  time  and  production  costs. 

The  disadvantages  of  the  programs  PRORUM  and  VORUM  are  that 
they  are  proprietary  and  were  specifically  developed  for  the  RUMX  2000 
machine  and  are  not  applicable  to  general  purpose  radial  forging  applications 
on  other  forging  machines. 

Metal  Flow  in  Radial  Forging 

The  metal  flow  of  solid  billets  in  radial  forging  has  been  analyzed  by 
Lange  (Ref.  1)  and  Paukert  (Ref.  21)  at  The  University  of  Stuttgart  using  FEM 
techniques.  Lange's  results  show  that  the  metal  flow  is  essentially  that  the 
material  underneath  the  forming  tools  will  be  displaced  as  a  rigid  block  into 
the  workpiece  while  the  regions  to  the  left  and  right  of  the  tool  will  be 
displaced  respectively  in  those  directions  with  no  corresponding  plastic 
deformation.  As  can  be  seen  below  in  Figure  12,  there  is  a  region  where  no 
axial  flow  of  material  occurs.  This  is  referred  to  as  the  neutral  plane.  Lange 
also  analyzed  the  conditions  which  allowed  maximum  plastic  deformation  in 
a  rectangular  cross-section  using  flat  tools  (Ref.  1).  To  obtain  maximum 
deformation  in  the  center  of  the  workpiece,  described  as  the  ratio  of  real  strain 
to  theoretical  strain,  Lange  found  that  this  occurred  when  the  ratio  of 
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indentation  to  original  billet  diameter  (bite  ratio)  was  equal  to  0.5.  To  ensure 
sufficient  deformation  in  the  workpiece  center,  a  relative  indentation, 
defined  as  the  ratio  of  2*Ah/H0,  of  at  least  0.16  was  required  (Ref.  1). 

Paukert  also  analyzed  the  radial  forging  process  using  visioplasticity 
and  FEM  for  a  flat  forging  tool  (Ref.  21).  Results  showed  that  the  most 
important  parameters  in  radial  forging  were:  ratio  of  tool  length  to  workpiece 
thickness,  workpiece  and  tool  geometry,  and  forging  sequence.  This  analysis 
needs  to  be  extended  to  the  case  of  tools  with  compound  angles. 

Isogawa,  Suzuki,  and  Uehara  more  recently  investigated  the  radial 
forging  of  Inconel  718  in  an  effort  to  study  the  nature  of  the  temperature 
profile  generated  and  the  resultant  microstructure  (Ref.  22).  As  part  of  their 
study,  plasticine  was  used  to  model  the  workpiece.  One  of  the  most 
interesting  points  brought  out  during  the  study  was  that  circumferential 
cracking  occurred  in  the  workpiece  at  a  critical  temperature.  Unfortunately 
the  authors  did  not  report  at  which  temperature  this  phenomena  was 
observed. 


Sb  =  Q9  1,2 

50mm/s 


Figure  12.  Velocity  Field  in  Radial  Forging  Showing  Effect  of  Relative 
Indentation  to  Core  Penetration  (Ref.  1). 
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Isogawa  et  al.  also  reported  that  the  shear  strain  increased  from  the  surface  to 
a  maximum  value  which  was  located  30%  below  the  surface  of  the  billet. 
After  reaching  a  maximum,  the  shear  strain  rapidly  dropped  to  a  negligible 
amount  with  increasing  distance  from  the  workpiece  surface.  This 
corresponded  to  the  region  at  which  the  authors  reported  that  circumferential 
cracking  had  occurred  during  actual  forging  trials.  The  cracking  effect  was 
attributed  to  the  radial  and  longitudinal  tensile  stresses  that  developed  in  the 
areas  that  were  between  those  in  contact  with  the  hammers. 

Cold  and  Warm  Forging 

As  the  trend  towards  net  shape  manufacturing  continues,  cold  and 
warm  forging  processes  offer  several  important  advantages:  reduced  energy 
consumption  for  billet  heating,  improved  surface  finish,  and  reduced  scale 
formation.  The  disadvantage  is  that  more  robust  forging  machines  are 
required  due  to  the  higher  flow  stresses  encountered.  The  primary  use  of  cold 
forging  in  radial  forging  is  limited  to  small  part  production  and  thin  wall 
tubes  due  to  the  higher  stresses  generated  during  deformation. 

Defect  Consolidation 

This  issue  is  of  considerable  importance  in  the  production  of  wrought 
billets  from  cast  ingots  in  primary  working  processes.  Defects  are  formed  from 
uneven  solidification  and  shrinkage  in  the  casting.  These  defects  seriously 
impair  the  soundness  and  mechanical  properties  of  the  final  forged  product. 
A  large  number  of  studies  have  been  performed  to  study  defect  consolidation 
in  open  die  forging  and  have  been  reviewed  in  an  earlier  ERC/NSM  report 
(Ref.  23).  However  this  issue  has  not  been  studied  as  extensively  for  the  radial 
forging  process.  It  would  be  of  practical  value  to  determine  how  the 
mechanism  of  void  closure  is  affected  by  four  radial  dies  as  compared  to  two 
dies.  Moreover  this  knowledge  would  be  of  use  in  optimizing  the  amount  of 
deformation  required  to  close  voids  which  are  present  in  the  billet.  Closing  of 
voids  in  open  die  forging  was  studied  by  Im  and  Dudra  (Ref.  23)  using  both 
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numerical  simulation  and  plasticine  modeling  techniques  for  plane  strain 
conditions  in  open  die  forging  of  massive  workpieces.  If  the  assumption  that 
the  effect  of  hydrostatic  stress  and  effective  strain  are  significant  as  claimed  by 
Im  and  Dudra,  their  approach  could  easily  be  extended  to  radial  forging.  Due 
to  the  fact  that  four  hammers  are  being  used  rather  than  two,  the  internal 
distribution  of  compressive  stresses  should  be  larger  and  hence  more  effective 
in  closing  internal  voids  and  porosity.  Although  this  is  intuitively  obvious,  it 
has  not  yet  been  formally  studied. 

Empirical  results  from  industry  have  shown  that  a  reduction  of  2.2:1 
will  consolidate  centerline  porosity  to  an  extent  that  it  can  no  longer  be 
detected  by  ultrasonic  inspection  in  wrought  bars  (Ref.  24).  An  increase  of  the 
reduction  to  3:1  has  been  found  to  achieve  absolute  core  density  and  a  fine 
grained  structure. 
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USE  OF  RFORGE  PROGRAM 


As  part  of  the  first  phase  of  this  work,  RFORGE  computer  program  was 
used  to  simulate  the  radial  forging  of  gun  tubes  at  The  Watervliet  Arsenal. 
RFORGE  is  presently  running  on  the  VAX  11/750  mini-computer  at  the  ERC 
for  Net  Shape  Manufacturing  at  The  Ohio  State  University.  A  user's  manual 
is  included  in  Appendix  A,  hence  only  an  overview  of  the  program  use  and 
capabilities  will  be  presented  here. 

Three  program  files:  INPUT.DAT,  OUTPUT.DAT,  and  RFORGE.FOR, 
are  required  to  run  the  simulation  routine.  All  input,  execution  and  process 
related  variables  used  in  the  simulation  are  entered  into  INPUT.DAT  and  are 
used  in  the  process  calculation  routines  contained  in  the  FORTRAN  77 
source  file:  RFORGE.FOR.  Upon  completion  of  program  execution,  all 
simulation  results  are  written  to  and  may  be  accessed  from  the  data  file: 
OUTPUT.DAT.  The  amount  and  detail  of  the  output  may  be  controlled  by 
inserting  the  appropriate  variable  value  in  the  input  data  file. 

RFORGE  is  essentially  a  fixed  configuration  FEM  model  for  forging  of 
hollow  tubes.  In  its  present  configuration,  RFORGE  is  limited  to  simulating 
forging  of  AISI  4337  tubular  workpieces  over  a  mandrel  except  in  the  case 
where  the  only  item  of  interest  is  the  forging  stresses  and  loads  for  which 
solid,  round  workpieces  may  be  used.  To  perform  the  simulation,  the  user 
defines  a  simplified  grid  pattern  for  the  workpiece  which  determines  the 
element  sizes  after  deformation  is  completed.  Deformation  is  assumed  to  take 
place  at  the  end  of  the  stroke.  As  no  plastic  deformation  is  taking  place  in  the 
sinking  zone,  only  the  forging  and  sizing  zones  are  considered  during  the 
analysis.  Workpiece  temperatures  are  calculated  by  determining  the  amount 
of  deformation  under  the  dies  assuming  no  heat  transfer  is  taking  place  at 
that  instant.  After  deformation  is  completed,  heat  transfer  is  then  assumed  to 
take  place  and  the  resultant  temperature  for  each  element  is  computed.  Stress 
and  load  calculations  are  performed  using  the  slab  method.  The  location  of 
the  neutral  plane  is  calculated  by  numerically  searching  each  of  the  three 
regions  for  the  point  in  the  velocity  field  where  the  material  has  zero 
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longitudinal  velocity.  Additionally,  this  region  also  corresponds  to  the 
maximum  radial  pressure. 

The  RFORGE  analysis  of  radial  forging  is  similar  to  that  of  extrusion. 
However  several  important  differences  exist  between  the  physical  processes 
and  must  be  taken  into  account  when  using  the  output  from  RFORGE.  As 
mentioned  earlier,  RFORGE  calculates  the  heat  added  to  the  workpiece  as  it 
occurred  at  the  end  of  the  stroke.  In  reality  the  ballet  is  deformed  to  its  final 
shape  in  a  series  of  increments  during  which  heat  generation  and  transfer  are 
taking  place.  As  such  the  heat  calculations  and  transfer  gradients  in  the  billet 
will  not  be  equal  for  each  approach.  RFORGE  in  this  case  would  be  an  upper 
bound  method  which  would  present  a  higher  temperature  distribution  than 
is  actually  present. 
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FINITE  ELEMENT  MODELING  OF  RADIAL  FORGING 


At  the  present  time  there  are  two  computer  programs  available  at  the 
ERC/NSM  which  may  be  used  to  analyze  the  radial  forging  process.  The  first 
program  is  a  dedicated  computer  program,  RFORGE,  which  was  developed  at 
Battelle  Columbus  Labs  under  contract  with  Watervliet  Arsenal.  The  second 
program  is  the  general  purpose,  FEM  based,  metal  forming  simulation 
package,  DEFORM,  which  was  also  developed  at  Battelle. 

For  the  following  simulations,  the  FEM  package,  DEFORM,  was  chosen 
over  RFORGE.  The  reasons  for  selecting  DEFORM  were: 

1.  Greater  flexibility  and  accuracy  to  model  a  wide  range  of  process 
conditions. 

2.  Ability  to  account  for  heat  loss  effects  not  directly  related  to  the  plastic 
deformation  process. 

3.  Capability  of  modeling  hollow  or  solid  workpieces. 

4.  Superior  post-processing  and  graphical  display  of  simulation  results. 

Initially  the  radial  forging  process  will  be  modeled  as  an  isothermal 
process  in  order  to  gain  increased  experience  in  the  use  of  DEFORM  for 
modeling  the  radial  forging  process.  The  isothermal  assumption  is  valid  due 
to  the  fact  that  in  the  forging  zone  that  heat  lost  to  the  environment  is 
replaced  by  the  heat  generated  from  the  dies  impacting  the  workpiece  rapidly. 
The  workpiece  temperature  may  in  fact  be  controlled  during  the  process  by 
varying  the  feed  of  the  workpiece  (Ref.  12).  However,  the  isothermal 
assumption  may  not  be  entirely  true  as  the  workpiece  temperature  may 
actually  increase,  and  hence  lower  the  flow  stress,  if  the  processing  conditions 
are  favorable.  However,  radial  forging  of  gun  steel  is  normally  done  in  the 
hot  regime  where  moderate  increases  in  temperature  are  not  expected  to 
significantly  impact  material  flow  stress  except  in  the  forging  of  titanium  and 
nickel  based  alloys. 
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Later  in  the  project,  the  process  will  be  modeled  non-isothermally  to 
include  the  effect  of  heat  loss  incurred  during  transport  of  the  billet  before 
and  during  the  deformation  process  to  study  the  effects  on  the  material  flow 
stress.  Additionally,  a  more  realistic  velocity  profile  will  be  incorporated  to 
account  for  the  characteristics  of  die  travel  similar  to  that  for  a  mechanical 
press.  As  the  software  becomes  available,  it  is  also  planned  to  model  the 
process  using  three-dimensional  FEM. 

The  results  presented  in  this  report  consist  of  two-dimensional, 
axisymmetric,  isothermal  simulations  for  AISI  4337  steel  using  concave  dies. 
Axisymmetry  in  the  model  implies  that  metal  flow  was  restricted  to  the  axial 
and  radial  directions.  Use  of  an  axisymmetric  model  minimized  the 
modeling  complexity  and  simulation  time.  Inherent  in  the  axisymmetric 
assumption  is  that  the  dies  completely  encircled  the  outer  surface  of  the 
workpiece  during  deformation.  As  such,  no  flash  was  extruded  between  the 
dies  in  the  model.  Additionally,  this  assumption  eliminated  the  need  to 
consider  the  workpiece  rotation  between  strokes  and  only  the  axial  feed  was 
considered.  To  accurately  model  the  process  using  flat  dies,  a  plane  strain 
model  or  three  dimensional  model  would  need  to  be  developed 
incorporating  the  workpiece  rotation.  Some  intermediate  and  final  results  of 
the  axisymmetric  simulation  are  shown  below  in  Figures  13-16.  The  process 
conditions  used  in  the  simulation  are  shown  below  in  Table  2. 
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Table  2.  Process  Conditions  Used  During  DEFORM  Simulation. 
Die  Velocity:  6.299  in/ sec.  (constant) 

Feed  Rate:  3.0  inches 
Material:  AISI  4337 

Billet  Temperature:  1652.0  degrees  F  (900  degrees  C) 

Heat  Transfer:  Isothermal 
Stroke  Increment:  0.415  inches 
Friction  Factor:  0.6  (Constant) 


Discussion  of  Results 

In  Figure  13  the  initial  mesh  and  the  rigid  die  used  throughout  this 
simulation  is  shown.  The  mesh  was  created  using  the  Automatic  Mesh 
Generator  (AMG)  and  system  mesh  optimizer  packages  in  DEFORM  which 
automatically  configures  an  optimum  mesh  based  on  the  workpiece 
dimensions.  For  this  simulation  a  22  x  4  node  mesh  was  generated.  To 
simulate  an  axisymmetric  workpiece,  the  symmetry  axis  was  oriented  along 
the  Y-axis  in  DEFORM.  This  enabled  the  simulation  time  to  be  optimized  as 
only  the  minimum  number  of  elements  are  included  in  the  mesh 
eliminating  redundant  computations.  For  simplicity,  a  flat  die  was  used  in 
the  simulation.  Boundary  conditions  imposed  on  the  workpiece  included  no 
metal  flow  across  the  centerline  of  the  billet  due  to  the  axisymmetric 
assumption  and  no  metal  flow  in  the  axial  direction  across  the  undeformed 
end  of  the  billet.  The  latter  restriction  was  based  on  the  assumption  that  in 
practice  that  the  'A'  chuckhead  would  be  gripping  the  workpiece  end  which 
will  prevent  metal  flow  at  the  end  of  the  workpiece. 

In  Figure  14  the  deformed  mesh  is  shown  with  the  corresponding  flow 
lines  generated  by  DEFORM.  The  neutral  plane  is  clearly  evident  in  the  center 
of  the  plastically  deformed  region.  Examining  the  flow  lines  at  the  neutral 
plane  show  that  all  metal  flow  is  directed  to  the  centerline  of  the  workpiece. 
Additionally  on  either  side  of  the  neutral  plane  the  flow  lines  show  that 
metal  flow  is  initially  directed  downward  and  then  flows  in  the  axial 
direction  as  expected  for  an  incompressible  material.  This  result  is  very  close 
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to  the  results  obtained  by  Paukert  (Ref.  21)  shown  in  Figure  12.  Also 
noticeable  is  the  effect  of  the  friction  factor  between  the  tool  and  workpiece 
surface.  It  may  be  seen  that  due  to  the  friction  conditions  specified,  some 
barreling  has  occurred. 

In  Figure  15  the  rigid  die  was  moved  upwards  by  3.0  inches  to  simulate 
the  axial  feed  of  the  workpiece.  Although  it  is  realized  that  in  practice  the  feed 
is  less  than  3.0  inches,  for  the  simulation  it  was  felt  that  a  smaller  feed  would 
not  yield  as  good  resolution  as  a  larger  feed  on  the  final  plots.  To  perform  this 
second  simulation,  the  deformed  mesh  from  Figure  14  and  duplicate 
deformation  and  boundary  conditions  were  imposed. 

Figure  16  shows  the  effective  strain  contours  which  resulted  from  the 
second  run.  The  contours  indicate  that  the  material  directly  underneath  the 
die  was  displaced  downward  as  a  rigid  block  of  material  causing  adjoining 
material  to  be  displaced  axially.  The  highest  strains  were  found  in  the  region 
lying  midway  between  the  die  and  the  workpiece  centerline.  This  region  was 
also  found  by  Isogawa  et  al.  (Ref.  22)  to  be  the  region  where  they  observed 
maximum  shear  strains  and  circumferential  cracking  in  their  study  of  radial 
forging  of  Inconel  718.  At  the  areas  directly  under  the  die  and  at  the 
centerline,  the  strains  were  found  to  be  approximately  50  percent  of  the 
maximum  effective  strain  values  demonstrating  that  deformation  was  to  the 
core.  An  interesting  result  is  shown  slightly  above  the  die  impression  where  a 
wavefront  pattern  of  strains  were  created  along  the  cross-section  due  to  the 
restricted  metal  flow  across  the  undeformed  workpiece  edge. 
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PRE-FORM  DESIGN  FOR  155  MM  GUN  TUBE 

One  problem  which  is  of  considerable  interest  in  radial  forging  is 
designing  preforms  which  will  yield  a  forged  product  which  is  near  net-shape 
and/or  has  consistent  mechanical  properties.  The  purpose  of  this  section  is  to 
study  the  experimental  data  which  was  provided  to  the  ERC  at  Ohio  State 
University  to  analyze  for  significant  factors  and  to  identify  an  optimum 
preform  for  the  155  mm  gun  tube  based  on  uniform  deformation  or  strain. 

The  focus  of  this  project  resulted  from  a  study  at  Benet  Laboratory 
which  evaluated  the  mechanical  properties  of  the  155  mm  gun  tube  for 
consistency  after  forging  and  heat  treatment.  The  properties  considered  in  the 
evaluation  were:  yield  strength,  tensile  strength,  percent  reduction  of  area, 
percent  elongation  and  Charpy  impact  values.  The  data  values  were  taken 
from  a  study  conducted  at  Watervliet  Arsenal  and  are  shown  below  in  Table 
3. 


Table  3.  Experimental  Data  From  155  mm  Gun  Tube  Forging. 


Total 

Forging 

Yield 

Tensile 

Charpy 

Reduction 

Reduction 

Strength 

Strength 

Impact 

(ksi) 

(ksi) 

(ft.-lbs.) 

12,  6  o'clock 

12,  6  o'clock 

3,  9  o’clock 

4.24:1 

1:1 

167.1, 164.4 

184.8, 183.0 

24.0, 23.0 

1.5:1 

164.4, 162.9 

183.0, 179.7 

21.0,24.0 

2:1 

165.6, 165.6 

183.0, 183.0 

23.0, 22.0 

10.73:1 

2.5:1 

165.3, 165.0 

183.0, 182.4 

23.0, 22.0 

12.91:1 

3:1 

165.0, 164.4 

182.4, 181.8 

22.0, 22.0 

15.10:1 

3.5:1 

165.0, 164.1 

182.7, 181.5 

22.0, 22.0 

17.31:1 

4:1 

164.4, 164.1 

181.2, 181.2 

21.0,23.0 

19.52:1 

4.5:1 

163.2, 163.2 

181.5, 181.2 

20.0, 21.0 

21.75:1 

5:1 

164.4, 164.4 

182.4, 182.4 

19.0, 20.0 
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The  most  critical  factors  as  identified  by  Engineers  at  Benet  Laboratory 
were:  yield  strength,  tensile  strength,  and  Charpy  Impact  Values.  Plots  of 
these  data  are  shown  below  in  Figures  17,  18,  and  19.  Examining  Figures  17, 
18,  and  19  shows  that  each  of  the  mechanical  properties  display  small 
amounts  of  scatter  when  plotted  against  forging  reduction.  It  is  also  apparent 
that  a  trend  exists  between  the  mechanical  properties  and  forging  reduction. 
To  study  the  exact  relationship  and  nature  of  the  trend  would  require  taking 
additional  data  points.  From  Table  3  the  range  of  the  data  was  found  to  be  3.9 
ksi,  3.6  ksi,  and  5.0  ft.-lbs.  for  the  yield  strength,  tensile  strength  and  Charpy 
Impact  respectively. 

To  determine  whether  the  mechanical  properties  in  Table  3  are  affected 
to  a  significant  extent  by  forging  reduction  and  heat  treatment,  it  was  decided 
to  test  the  data  using  Analysis  of  Variance  Techniques  (referred  to  as 
ANOVA).  ANOVA  methods  are  powerful  statistical  methods  and  procedures 
used  to  analyze  experimental  data  (Ref.  26).  For  the  ANOVA  analysis 
presented,  a  single  factor  was  tested  to  determine  whether  or  not  it  was 
affected  by  forging  reduction  and  heat  treatment  or  if  the  variation  was  due  to 
inherent  randomness  in  the  data.  Additionally,  the  difference  in  location  was 
not  considered  as  a  variable  but  rather  as  a  replication.  The  specific  ANOVA 
method  employed  was  the  single  factor  ANOVA  test.  Essentially,  the  single 
factor  ANOVA  is  testing  two  different  hypotheses  using  the  well  known  F- 
Test  to  determine  whether  or  not  there  are  differences  in  the  true  averages  of 
the  different  populations  due  to  the  treatment  (Ref.  26).  As  the  relevant 
theory  and  equations  are  presented  in  many  statistical  texts,  only  the  final 
calculated  results  will  be  presented  here.  The  resulting  ANOVA  tables  are 
shown  below  in  Tables  4,  5,  and  6.  For  the  F-test,  a  significance  level  of  0.05 
was  used  with  eight  degrees  of  freedom  (referred  to  as  d.o.f.)  for  the 
numerator  and  nine  d.o.f.  for  the  denominator.  The  critical  value,  F (0.05,8,9), 
obtained  from  a  table  of  F-test  values,  for  rejecting  the  null  hypothesis  was 
3.23  for  all  three  tests.  The  null  hypothesis  for  the  experiment  states  that:  the 
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Figure  17.  Yield  Strength  Plots  for  AISI 4337  Steel  Data. 
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Figure  19.  Charpy  Impact  Values  for  AISI 4337  Steel. 


means  of  the  data  at  each  forging  reduction  are  equal.  Equivalently,  it  may  be 
stated  that  there  is  no  difference  in  a  particular  mechanical  property  as  forging 
reduction  is  changed  in  the  test  under  consideration.  The  alternative 
hypothesis  states  that  the  means  of  the  data  are  not  equal,  hence  the  given 
mechanical  property  is  affected  by  forging  reduction. 


Table  4.  ANOVA  Table  for  Yield  Strength  Test. 


Source  of 

D.O.F. 

Sum  of 

Mean 

F 

Variation 

Squares 

Square 

Reduction 

8 

11.320 

1.415 

2.339 

Error 

9 

5.445 

0.605 

Total 

17 

16.765 

Table  5.  ANOVA  Table  for  Tensile  Strength  Test. 


Source  of 

D.O.F. 

Sum  of 

Mean 

F 

Variation 

Squares 

Square 

Reduction 

8 

12.550 

1.569 

1.724 

Error 

9 

8.190 

0.910 

Total 

17 

20.740 

Table  6.  ANOVA  Table  for  Charpy  Impact  Test. 


Source  of 

D.O.F. 

Sum  of 

Mean 

F 

Variation 

Squares 

Square 

Reduction 

8 

22.778 

2.847 

2.847 

Error 

9 

9.000 

1.000 

Total 

17 

31.778 
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As  can  be  seen  in  Tables  4,  5,  and  6,  none  of  the  F-test  values  were 
found  to  be  greater  than  F(0.05,8,9).  Hence  the  null  hypothesis  could  not  be 
rejected  in  any  of  the  cases  under  consideration  using  the  given  data.  That  is 
to  say  that  yield  strength,  tensile  strength,  and  Charpy  Impact  Values  are  not 
greatly  affected  by  the  combined  forging  and  heat  treatment  process  being 
used.  This  does  not  totally  eliminate  the  possibility  that  forging  reduction 
does  affect  mechanical  properties  due  to  the  large  error  terms  found  during 
the  analysis.  If  the  error  terms  could  be  reduced  then  the  analysis  would  show 
that  there  is  a  significant  effect  between  forging  reduction  and  the  mechanical 
properties.  A  large  error  term  indicates  that  a  large  amount  of  random  error  is 
contained  in  the  data.  However,  based  on  the  results  of  the  ANOVA  tests  in 
Tables  4,  5,  and  6,  it  is  not  possible  to  design  an  optimum  preform  from  the 
given  data.  Instead  the  problem  requires  additional  study  to  eliminate  the 
confounding  effects  of  one  process  upon  the  other.  To  accurately  determine 
whether  the  radial  forging  process  affects  the  mechanical  properties  of  AISI 
4337  under  the  present  forging  conditions,  a  separate  study  is  required.  A 
recommended  approach  is  to  forge  at  temperature  specimens  of  AISI  4337  at 
differing  amounts  of  strain,  air  cool,  and  perform  the  required  mechanical 
tests.  The  next  step  should  be  to  separately  investigate  the  effects  of  the  heat 
treatment  process  on  mechanical  properties  as  material  volume  or  forging 
reduction  is  varied  using  the  same  forging  conditions  in  use  at  Watervliet.  It 
is  recommended  that  these  issues  be  further  studied  prior  to  changing  the 
present  preform  design.  This  is  an  important  issue  since  the  forging  process  is 
being  carried  out  in  the  warm  working  range  and  the  relative  effects  of  work 
hardening  and  recovery  on  the  resultant  microstructure  is  not  yet  known. 

Preform  Design  for  Desired  Product  Properties 

In  the  thermomechanical  processing  of  alloys  there  are  two  basic 
philosophies  that  may  be  employed  to  yield  a  desired  microstructure.  The  first 
approach  is  to  produce  the  desired  shape  by  any  means  possible  and  then  use 
heat  treatment  to  produce  the  desired  microstructural  properties  (Ref.  27). 
The  second  approach  is  to  tailor  the  mechanical  processing  at  achieving  the 
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preselected  microstructure  and  abiding  by  the  processing  limitations  imposed 
on  temperature  and  percent  reduction  (Ref.  27).  In  this  study  the  primary 
objective  is  to  further  the  goal  of  achieving  a  desired  microstructure  during 
the  radial  forging  process.  Hence  the  latter  approach  will  be  employed. 

This  study  is  an  effort  to  ensure  that  uniform  mechanical  properties 
are  produced  in  forged  gun  tubes  through  application  of  uniform  plastic 
deformation  during  the  forging  process.  The  primary  focus  is  on  Charpy 
Impact  Data  which  indicates  the  relative  toughness  of  the  workpiece.  In 
general,  the  toughness  of  a  workpiece  decreases  with  increasing  amounts  of 
cold  work.  However,  production  of  gun  tubes  is  normally  carried  out  in  the 
warm  or  hot  working  regions  where  the  workpiece  undergoes  recovery  and 
recrystallization  during  deformation.  For  this  study  a  preliminary 
investigation  of  flow  stress  curves  at  warm  working  temperatures  from  a 
prior  study  by  Battelle  Columbus  Labs  (Ref.  25)  showed  a  moderate  increase  in 
flow  stress  with  increased  deformation.  It  is  well  known  that  toughness  is  a 
function  of  grain  size.  For  very  small  grain  sizes,  toughness  increases  with 
grain  size  to  a  maximum  value  after  which  toughness  decreases  with  further 
increases  in  grain  size.  Hence,  toughness  is  normally  gained  at  the  expense  of 
strength  under  cold  working  conditions  except  when  it  is  achieved  through 
grain  refinement  (Ref.  28). 

Grain  Growth  and  Recrystallization  Kinetics 

Before  describing  the  effects  of  deformation  on  the  recrystallization 
process  it  is  first  helpful  to  describe  recovery  and  recrystallization.  Recovery  is 
a  process  which  acts  to  reduce  the  strain  energy  of  the  deformed  material  to  a 
more  thermodynamically  stable  state  by  eliminating  both  vacancies  and 
dislocations  without  any  subsequent  rearrangement  of  the  existing 
microstructure.  As  such,  the  original  grain  size  is  retained  after  recovery 
occurs.  Another  point  that  must  be  noted  is  that  recovery  occurs  at  a  lower 
temperature  than  does  recrystallization.  Recrystallization  also  acts  to  reduce 
the  strain  energy  imparted  by  deformation  but  accomplishes  this  through 
nucleation  and  growth  of  new,  strain-free  grains.  This  process  occurs  at  a 
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higher  temperature  than  does  recovery.  Before  recrystaiiization  can  take 
place,  a  minimum  amount  of  deformation  or  percent  reduction  must  be 
imparted  to  the  material. 

The  two  processes  may  be  considered  to  be  in  competition  with  one 
another  in  determining  the  final  grain  size.  Once  recovery  or  recrystallization 
has  taken  place,  then  it  is  not  possible  for  the  other  process  to  occur.  Both 
processes  are  dependent  on  the  rate  of  heating  as  well  as  temperature  (Ref. 
29).  If  a  slow  rate  of  heating  is  used  while  heating  to  the  recrystallization 
temperature  then  more  recovery  is  able  to  occur  as  additional  time  is  spent  in 
the  recovery  region  and  more  material  will  be  reclaimed  by  this  process.  The 
reverse  also  applies  to  recrystallization. 

It  is  well  known  that  deformation  affects  the  recrystallized  grain  size  of 
most  steels.  Under  cold  working  conditions  the  recrystallized  grain  size  will 
decrease  as  percent  strain  increases  during  annealing.  Therefore,  it  is  clear 
that  forging  reduction  is  an  important  factor  in  controlling  the  resultant  grain 
size.  The  initial  grain  size  of  the  preform  also  influences  the  size  of  the 
recrystallized  grain  (Refs.  28,29).  This  may  be  seen  by  noting  that  new, 
recrystallized  grains  are  nucleated  at  primarily  at  grain  boundaries,  triple 
junction  points  of  boundaries,  precipitates,  and  other  points  where  regions  of 
inhomogeneous  plastic  deformation  exist  to  create  high  concentrations  of 
strain  energy  which  drive  grain  growth.  In  particular  if  a  smaller  grain  size  is 
present  in  the  preform,  more  nucleation  points  will  be  available  for 
recrystallization.  To  effectively  control  the  final  grain  size  of  the  gun  tube  and 
to  optimize  the  forging  process,  it  is  necessary  to  determine  at  what 
temperatures  that  recovery  and  recrystallization  occur  for  AISI  4337. 

Effect  of  Deformation  Parameters  on  Recrystallization 

The  importance  of  strain  and  strain  rate  on  recrystallization  must  not 
be  ignored  as  they  directly  impact  the  final  grain  size.  The  importance  of 
strain  rate  may  be  described  as  follows.  A  larger  strain  rate  enables  higher 
processing  or  production  rate.  Heat  loss  to  the  dies  and  workpiece  chilling  is 
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also  minimized  (Ref.  28).  The  result  of  heat  loss  or  chilling  will  obviously 
lower  the  billet  temperature  and  cause  a  corresponding  increase  in  the  flow 
stress.  In  the  radial  forging  process  this  would  occur  primarily  at  the  inner 
and  outer  diameters  of  the  tube  at  the  tube/mandrel  and  tube/die  interfaces. 
Moreover,  this  could  lead  to  a  non-uniform  distribution  of  deformation 
across  the  tube  wall  which  could  ultimately  affect  the  recrystallized  grain 
structure.  It  is  not  anticipated  that  workpiece  chilling  due  to  die  contact  is 
significant  due  to  the  relatively  short  contact  time.  Heat  loss  by  convection 
and  radiation  to  the  environment  and  to  the  water-cooled  mandrel  will  be 
more  significant.  However,  the  heat  losses  may  be  compensated  through 
proper  selection  of  the  stroking  rate  and  feed  rate.  This  is  also  favorable  in 
that  increased  rate  of  heating  minimizes  recovery  during  deformation. 
However,  the  penalty  for  increased  strain  rate  is  that  flow  stress  will  also 
increase.  This  may  be  seen  in  the  following  well  known  power  law  relation: 

(10)  a  =  Ce 


The  effect  of  strain  may  also  influence  the  temperature  of  the 
workpiece.  Intuitively  as  the  amount  of  deformation  is  increased,  greater 
amounts  of  work  must  be  supplied  to  the  billet.  In  plastic  deformation 
processes,  95%  of  the  added  work  will  show  up  as  heat  in  the  form  of 
temperature  increase.  Mathematically  the  change  in  temperature  may  be 
described  as: 

(11)  To  -  aA£P  (Ref.  7) 

0  Jcp 

Hence  larger  strains  will  also  tend  to  promote  increases  in  temperature  as 
does  strain  rate. 

Effect  of  Grain  Refiners  on  Recrystallization 

Another  possibility  which  exists  for  controlling  the  grain  size  in  the 
forging  is  to  investigate  the  introduction  of  stable  second  phase  particles  or 
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grain  refiners  into  the  steel  matrix.  Second  phase  particles  are  commonly 
used  to  control  grain  size  in  aluminum  killed  steels  and  aluminum  alloys. 
The  grain  refiners  are  most  effective  when  they  are  dispersed  in  a  fine,  even 
distribution  throughout  the  matrix  where  they  act  as  a  drag  against  expanding 
grain  boundaries.  In  order  for  the  grain  boundaries  to  continue  expanding, 
additional  energy  must  be  available  to  the  grain  boundary  (Refs.  28,29).  The 
effect  of  the  second  phase  particle  size  on  final  recrystallized  grain  diameter 
may  be  seen  in  the  Zener-Holoman  equation  below.  Caution  must  be 
exercised  when  the  material  is  heated  to  avoid  the  reversion  temperature  of 
the  particles  at  which  they  will  dissolve  back  into  solution  and  permit 
discontinuous  grain  growth  to  take  place. 

(12)  D  (Ref.  28) 

3v7 

Preform  Development 

The  present  preform  shape  being  used  at  Watervliet  is  shown  below  in 
Figure  20.  The  new  preform  design  will  be  developed  using  the  concept  of 
constant  volume  during  plastic  deformation.  The  concept  of  constant  volume 
in  plastic  deformation  processes  is  based  on  the  fact  that  when  all  strain 
components  are  added  together  the  result  is  zero  which  indicates  that  no 
change  in  workpiece  volume  has  occurred.  This  is  written  mathematically  in 
Cartesian  coordinates  as: 

(13)  e„+e„  +  eu=0 

Here  £»  is  defined  as  true  or  logarithmic  normal  strain  as  in  equation  14: 

.  final  dimension 

(14)  e„  =  In  —  — - — 

initial  dimension 

Although  the  subscripts  xx,  yy,  zz  are  commonly  used  in  many  forming 
analyses,  the  polar  coordinate  system  is  more  commonly  employed  in  studies 
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involving  prismatic  or  rotational  workpieces  involving  circular  and/or 
tubular  cross-sections.  The  notation  for  polar  coordinates  is  as  follows:  £n,  £qq, 
£zz  corresponding  to  the  radial,  tangential,  and  axial  directions  as  shown 
below  in  Figure  21.  The  strains  are  calculated  as  follows  (Ref.  30): 

(15) 


(16) 


e_  =  In 

rr  1 

=  In 
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Figure  20.  Experimental  Preform  Design  (Ref.  Sketch  7-89). 


Figure  21.  Polar  Coordinate  Notation. 
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(17)  £22  =  -(Err  +  e96^  ^rom  constant  volume 
alternatively  Ea  may  be  expressed  as  : 

(18)  £„  =  hA 

z* 

To  facilitate  the  analysis,  the  finish  forged  gun  tube  was  divided  into 
three  geometrically  distinct  regions  as  shown  below  in  Figure  22.  It  was 
planned  to  calculate  the  strains  from  the  muzzle  end  of  the  gun  tube  as  the 
data  indicated  that  the  highest  Charpy  Impact  values  were  found  at  the 
muzzle.  This  implied  that  as  for  the  as-received  preform,  that  minimum 
deformation  and  the  subsequent  heat  treatment  yielded  the  highest 
toughness.  However  upon  analyzing  the  strains  at  the  muzzle  end  it  was 
found  that  due  to  the  small  amount  of  deformation  applied,  the  wall 
thickness  had  increased  (i.e.  positive  value  of  £„).  This  condition  is  often 

observed  in  tube  sinking  for  low  values  of  friction  and  reduction  (Ref.  30). 
However  as  these  conditions  would  be  applicable  only  during  deformation  at 
the  muzzle  end ,  applying  the  same  conditions  throughout  the  workpiece 
would  not  yield  a  satisfactory  design  for  the  preform.  Therefore,  it  was 
decided  to  use  the  reduction  at  the  breech  end  to  calculate  the  strains  for  the 
preform  as  this  was  felt  to  be  more  representative  of  the  overall  forging 
conditions.  It  should  be  noted  that  a  different  preform  could  be  designed  to 
optimize  the  microstructure  for  recrystallization  by  simply  replacing  the 
values  of  £„,  Eee,  and  with  the  desired  values  of  reduction  and  repeating 

the  calculations  shown  below.  It  was  assumed  during  the  analysis  that 
sufficient  deformation  would  be  available  during  forging  to  drive 
recrystallization.  However  this  should  be  verified  with  the  material  vendor 
as  this  information  was  not  available. 
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Figure  22.  Definition  of  Geometric  Regions. 

Calculating  the  strain  at  the  breech  end  the  calculations  are  shown  below: 
6.225  -  2.660 

(19)  £„  »  In  7  5QQ .  3  250  *  ' 0 176 


6.225  +  2.660 

(20)  £e0  =  In  7  500  +  3  250  =  "°*191 


(21)  £ ^  *  -(-0.178  -  0.191)  =  0.367 


Now  applying  this  to  the  preform  to  determine  the  axial  length,  Zj  of  the 
preform: 


<22> 


exp(0.367) 


Using  zf  =  82.5  inches  from  Figure  20,  z.  =  57.157  inches. 

To  perform  a  similar  analysis  for  the  muzzle  end,  it  must  be  noted  that 
to  solve  for  the  preform  dimensions  that  the  equation  is  a  function  of  two 
variables:  Rq  and  Rj.  This  is  solved  by  noting  that  the  internal  bore  for  the 
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preform  will  be  the  same  throughout:  6.500  inches.  Using  a  similar  analysis  as 
in  equation  19  for  the  breech  end  of  the  preform  the  following  values  are 
obtained: 


(23) 


4.025  -  2.660 
n  R0-  3.250 


-0.176 


solving  for  RQ/  Rq  =  4.880  inches.  Similarly  for  the  preform  length: 

82.500 

(24)  In - =0.367 

z. 

l 

which  leads  to  the  result  that  z.  =  64.431  inches. 

i 

For  Region  n  the  preform  dimensions  can  be  calculated  by  noting  the 
fact  that  the  volume  of  the  preform  and  the  forging  are  equal  and  that  the 
inner  and  outer  radii  of  the  preform  at  the  entrance  and  exit  of  Region  II  are 
equal  to  those  of  Region  I  and  III  respectively.  Hence  the  problem  reduces  to 
that  of  equating  the  volumes  of  the  preform  and  the  workpiece  and  solving 
for  Zj.  The  equation  for  the  volume  of  a  frustum  is: 


(25) 


7t[r2  +  rR  +  R2]h 
Volume  =  - - - 


where  r,  R,  and  h  are  shown  in  Figure  23.  To  account  for  the  internal  bore  of 
the  workpiece,  the  volume  of  a  right  circular  cylinder  is  subtracted: 

(26)  Volume  =  Jtr2l 

where  r  is  the  radius  of  the  cylinder  and  1  is  the  length. 

Writing  the  equation,  volume  preform11*  volume  forging11: 
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Figure  23.  Notation  for  Calculation  of  Frustum  Volume. 

(27)  (7.5002  +  4.880x7.500  +  4.8802)  z.  ^  -  tc(3.250)2Zi  = 

(6.22 52  +  4.025x6.225  +  4.0252)3(£  -  ji(2.660)230 

3 

Zj  is  found  to  be:  20.752  inches.  The  resulting  preform  design  is  shown  below 
in  Figure  24. 
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Figure  24.  New  Preform  Design  for  Watervliet  Arsenal. 
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CONCLUSIONS 


In  comparison  with  other  forging  processes,  the  state  of  radial  forging 
process  analysis  has  not  been  highly  developed.  Of  the  forming  models  which 
have  been  presented,  most  rely  on  the  approximation  that  the  process  may  be 
modeled  as  a  variation  of  extrusion.  While  this  approximation  greatly 
facilitates  the  mathematical  modeling  of  radial  forging,  it  does  not  accurately 
reflect  the  metal  flow  and  friction  conditions  which  have  different  signs  on 
either  side  of  the  neutral  plane  as  in  upsetting.  With  the  advent  of  FEM 
programs  such  as  ALPID  and  DEFORM,  radial  forging  process  modeling  has 
been  greatly  facilitated  and  all  conditions  are  easily  simulated,  including 
forging  of  hollow  and  solid  workpieces.  However,  a  true  simulation  of  the 
physical  process  can  not  be  made  until  new  three-dimensional  FEM  software 
becomes  available  and  accurate  flow  stress  data  for  metals  at  warm  and  hot 
working  temperatures  is  generated. 

1.  In  this  study  RFORGE  and  DEFORM  were  used  to  perform  radial  forging 
simulations.  DEFORM  was  found  to  be  superior  to  RFORGE  for  simulating 
a  wide  variety  of  process  conditions  but  neither  was  found  to  be  useful 

in  preform  design. 

2.  Deformation  conditions  during  radial  forging  have  a  direct  effect  on  the 
resulting  microstructure  of  the  forged  tube.  However  due  to  probable 
confounding  effects  from  heat  treatment,  no  effect  was  noted  on  the 
mechanical  properties  with  varying  reduction. 

3.  The  effects  of  the  radial  forging  and  heat  treatment  processes  should  be 
studied  independently  to  ascertain  their  effect  on  the  microstructure  of 
the  gun  tube. 

4.  For  the  preform  design,  simple  plastic  strain  calculations  and  concepts 
were  used  to  develop  the  design  presented  in  Figure  24.  It  is 
recommended  that  the  true  strain  concept  be  adopted  for  describing 
the  forging  reduction. 
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APPENDIX  A:  RFORGE  USER’S  MANUAL 
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DESCRIPTION  OF  THE  COMPUTER  PROGRAM  RFCRGE 


During  Phase  I  and  Phasa  II  of  this  dsveloposnc  scudy  on  the 
"Analysis  and  Optimisation  of  tha  Radial-Forging  Process  for  Manufacturing 
Cun  Barrels",  three  sets  of  computer  programs  were  developed.  The  first 
set,  named  RADFLW,  determines  a  kinematically-admissible  upper-bound 
velocity  field,  strains,  strain  rates,  and  an  upper  bound  on  the  total 
forging  load  in  radial  forging  of  gun  barrels  over  a  stationary  mandrel. 

The  second  sec  of  programs  called  RFTEMP ,  uses  the  velocity  field  data 
generated  by  RADFLW  to  obtain  temperature  distributions  in  the  billet,  the 
mandrel,  the  die  and  the  product  in  radial  forging  of  tubes.  The  third  sec 
of  programs,  named  RFCOMP ,  determines  the  stresses  and  loads  in  radial  forg¬ 
ing  of  gun  barrels  or  cubes  using  the  compound  angles,  and  it  is  also  valid 
for  dies  with  single-entrance  angle. 

For  overall  optimization  of  the  radial-forging  process,  all  these 
three  sets  of  programs  are  grouped  together  under  the  name  RFORGE.  In  this 
modified  form,  any  desired  information,  such  as  upper-bound  velocity  field, 
strains,  strain  rates,  temperature  distributions  or  the  load  and  pressure  on 
the  dies,  can  be  generated  with  one  sec  of  input  data.  Thus,  with  the  use  of 
RFORGE,  generating  data  to  plot  performance  curves  is  relatively  simple,  and 
the  process  optimization  for  a  given  set  of  -nput  conditions  can  be  easily 
established  with  a  minimum  number  of  computer  runs. 

All  the  routines  in  RFORGE  have  been  coded  as  subroutines  or  sub- 
functions.  Bence,  they  can  be  used  in  a  simple  one-line  structure  or  in  an 
overlay  structure.  Since  the  memory  space  requirements  for  the  one-line 
structure  of  the  RFORGE  exceeded  the  capacity  of  the  IBM-360  computer  at  the 
tfatervliec  Arsenal,  an  overlay  structure,  as  depicted  in  Figure  B-l,  was 
developed.  In  its  present  overlay  form,  RFORGE  requires  19C88  (hexadecimal) 
or  81000  (decimal)  memory  locations.  Since  the  overlay  structure  is  established 
by  a  set  of  control  cards,  modification  of  the  structure  due  to  the  addition  or 
deletion  of  subroutines  requires  minimum  effort. 
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RFOROE 
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FIGURE  B- 1 .  OVEKI.AY  STRUCTURE  OF  TUE  SET  OF  COMPUTER  PKOCKAMS  RFOROE 


All  process  variables  are  '■cad  by  RfQRGE.  In  addition,  the  choice 
co  compute  one  or  more  of  (a)  upper-bound  velocity  field,  (b)  temperature 
distributions  in  the  billet,  the  mandrel,  the  die  and  the  product,  and  (c) 
load  and  pressure  on  the  dies,  are  read  in  as  input.  The  process  variables 
being  read  include: 

(a)  The  dimensions  of  the  preform  and  of  the  forged 
barrel 

(b)  The  dimensions  and  the  velocities  of  the  forging 
dies 

(c)  The  values  of  friction  at  die  material  and  mandrel- 
material  interfaces 

(d)  The  magnitudes  of  back  push,  front  pull,  and  axial 

feed 

(e)  Initial  temperature  and  type  of  the  metal  being 
forged. 

The  flow  stress  of  the  forged  material  for  given  values  of  strain,  strain  race, 
and  temperature  are  obtained  by  interpolation  from  the  tabulated  flow  stress  da 
Included  in  the  program. 

RFORGE,  which  coordinates  the  functions  of  all  ocher  major  routines, 
also  prints  out  the  data  read  in.  After  the  initialization  of  some  of  the  inpu 
values,  the  major  routines  RADFLW,  RETT*? ,  and  ?FC0M?  are  called  as  prescribed 
by  the  user  during  the  input  procedure. 

(*) 

Input  Preparation  for  RFORGS 

(The  quantity  la  parentheses,  following  the  description  of  each  input  variable, 
represents  an  example  Input) . 

(TRUE  or  FALSE)  input  to  indicate  whether 
velocity  field  Is  to  be  generated.  A  'T'  in  anyone  of 
,  the  first  10  columns  will  set  this  variable  co  TRUE. 

An  ’F*  or  ell  blanks  will  set  it  to  FALSE.  (T) 

(*)  See  Figure  B-2  for  variable  names  used  here. 


1.  First  Card,  Format  (3L10.110] 
Columns  1-10:  A^  \  Lc 


To  be  included 
only  when  "c" 
in  the  first  cord 
is  sel  to  TRUE 


PI CURB  U-2 ,  INPUT  DATA  STRUCTURE 


11-20:  j.  Logical  (T3.UZ  or  FALSE)  inout  Co  indicate 
whether  cha  alaboraca  temperature  distri¬ 
butions  ara  Co  ba  computed.  Sine*  velocity 
field  la  naadad  for  cha  computation  of 
caaparacura  distribution,  cha  first  variable 
ausc  also  ba  sac  co  crua,  whan  chls  variable 
is  assigned  crua.  ’.fhan  cha  first  variable  is 
sac  co  TRUE  and  cha  second  variable  is  set  to 
FALSE,  an  average  temperature  of  the  deforma¬ 
tion  zone,  heat  generated  by  metal  deformation, 
haac  absorbed  by  deforming  material,  heat  lost 
by  conduction,  convection  and  radiation  are 
calculated  by  approximate  heat  transfer  formiias  .  (T) 
21-30:  C,  Logical  (TRUE  or  FALSE)  input  to  indicate  whan 

load  and  pressure  estimations  are  to  be  performed.  (7) 
31-40 :  D,  An  integer  used  to  print  out  intermediate  results. 
When  set  to  zero,  only  the  final  results  will  be 
printed.  Input  data  and  final  results  will  be 
printed  when  set  to  1.  All  other  variables  that 
will  be  printed  for  a  higher  setting  of  this 
variable  are  shown  in  Table  3-2  attached  at  the 
end.  Maximum  value  for  IPRINT  is  9.  Faulty 
setting  of  10  or  more  will  automatically  be 
reduced  to  1.  (3) 


2.  Second  Data  Card,  Format  (3F10.3) 

Column*  1-10:  E,  Outside  diameter  of  the  tubular  preform,  mm  (560.) 

F,  Inside  diameter  of  the  tubular  preform,  cn  (313.) 

G,  Outside  diamscar  of  the  forged  product,  mm  (353.) 

H,  Inside  diameter  of  the  forged  product,  mm  (287.) 

I,  Length  of  the  die  land,  mm  (165.) 

J,  Radial  velocity  of  the  forging  cools,  mm/sec  C130.) 

K,  Axial  feed  of  the  preform,  m/stroke  (10.) 

L,  Entrance  angle  of  the  forging  cools,  degrees  (15.) 


11-20: 

21-30: 

31-40: 

1-50: 

1-60: 

:l-70: 

71-80: 
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Column  1-10:  M,  Back-pull  fore*  (negative,  if  back  push),  kN  (0.) 

11-20:  N,  Front-pull  fores  (negative,  tf  front  push) ,  kN  (0.) 
21-30:  0,  Friction  shsar  factor  between  tuba  and  dla  (0.6) 

31-40:  P,  Friction  shsar  factor  batvssn  tubs  and  aandrel  (0.5) 
41-50:  Q,  Initial  unifora  tsapsraturs  of  ths  billst ,  C  (1000.) 
51-60:  R,  (3)  angle  of  contact  for  a  round-facsd  tool,  degrees. 

(Sss  Figure  3-3,  optional  for  a  flat-facsd  tool)  ( 90 . % 
61-70:  S,  (;)  angls  of  cangency  of  a  flat-facsd  tool,  dsgrsss. 

(Sss  Figurs  3-4,  optional  for  a  round-facsd  tool)  (0.) 
71-30:  T,  A  cods  for  type  of  tool  (1) 

1  for  round-facsd  tool 

2  for  flat-facsd  cool. 

As  an  example,  ths  input  data,  as  princsd  by  ths  program  RFCRGE,  ars 
show  in  Figurs  3-5.  Ths  program  prints  chasa  input  data  so  that  input  can  bs 
chacksd  for  posalbla  srrors. 

RF0RGZ  rsads  ths  data  continuously  until  ths  and  of  fils  is  encountered. 
Hanes,  for  procsss  optimization,  any  nuabsr  of  data  sacs  can  bs  read  in  a  single 
run. 

The  following  data  card  should  not  bs  included  in  ths  input  data  stream, 
if  ths  first  variable  (A)  in  the  input  list  to  3J0RGE  is  sac  to  FALSE. 


4,  Fourth  Card,  Foraat  (3110) 

Columns  1-10:  0,  Number  of  radial  divisions  on  the  entrance  side  (15) 

11-20:  7,  Number  of  axial  divisions  (12) 

21-30:  V,  Number  of  radial  divisions  on  the  exit  slds  (5) 

As  an  example,  the  input  data,  as  princsd  by  the  program  RADFLV,  ars 
given  in  Figure  B-6.  The  sample  output  for  the  same  input  data  is  given  in 
Figure  8-7. 
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input  to  rfor 


.  variable  to  execute  velocity  program. r  ; 


QGICAL  VARIABLE  TQ  EXECUTE  TEMPERATURE  PROGRAM. F  m 


SI  CAL  VARIABLE  TO  EXECUTE  PPESSUPE  PPQGP  AH  .F  • 


NOEX  TO  POINT  INTERMEDIATE  RESULT 


OUTSIDE  OIAMETER  OF  THE  Tu9ULAR  PREFORM.mm  • 


.AR  PREFORM  .MM 


OUTSIDE  OIAMETER  OF  THE  FORGED  PRODUCT  .MM 


INSIOE  diameter  OF  TM£  FORGED  PRQOuCT.mm 


FIGURE  B-5.  SAMPLE  INPUT  DATA  AS  PRINTED  BY  RJORGZ 
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TNP»1T  TO  The  PROGRAM  PAOFLW 


number 

OF 

RAO  I AL 

UIVISIONS 

ON 

ENTRY  SIOE 

a 

15 

number 

OF 

total 

axial  oivisions 

a 

12 

NUMBER 

"of“ 

RAO I AL 

OIVISIONS 

ON 

EXIT  SIDE 

a 

5 

FIGURE  8-6.  SAMPLE  INPUT  DATA  AS  PRINTED  BY  RADFLtf 
PRIOR  TO  PERFORMING  COMPUTATIONS 


69 


_CUTO.jr  PROrt  THE  PROGRAM  SAQFLW 


LENGTH  OP  TM£  PURGING  ZC‘Ig,MH  a 
OISTANCE  op -iHt  NEtjTRA  l.  PLANE  PRQM t NT  RyTmN^1 
AVERAGE  STRAIN  In  Th£  PCRGING  ZOSE.m/h  a 
AVERAGE  ST(\AIN-SAT*  in  THc  PORGING  ZQN£,i/S  a 
AicRAGE  PLCW  S T F23S *  “M/SO.  M  '  5" 
TOTAL'  PAOIAL  PORGING  L0A3  PER  TQCL.N  a 
AVERAGE  RACIAL  PCRGI NG  PR£ SSUR £, “M/SG. N  a 
ENERGY  RATE  SJPBLIE3  9 Y  TQCLS.nh/S - T 


. ?36267E*Q3 
TMt2T7ToJ 
.  17;gG52*Q1 

•  13  5453£*Jl" 

T  i  iT93TSTqT 

•  45l932i*03* 
.2«.530l£*03' 

nwmrrrr 


G^(J)  »  AXIAL  COORDINATE  OF  A  C  =  IC  POINT.NM 

/T2  /  T  1  «  a  «  ^  •  •"  "  «**!•"  "T  *  _"“T_  • -  *> 


GR(I)  a  RACIAL  CCGPG  INATE  "oT' A’  GRlO--?*!  NTThTT 
..UZ(I,J)  a  AXIAL  CCHOQNENT  OF  VfLOC I7Y , in/ SEC 
U  R  ( I  »  J  )  a  RACIAL  CQ-FONcNT  CP  VLLCC  IT  Y  f  hh/ScC* 
_ SR(I,JI  a  cPfrr-rivE  SYRAJ  N- FA  T£»t/S  EC 


-I _ 

J 

GZtJ) 

_  G  S  ( I ) _ 

UZCI.J) 

URII.J) 

1 

1 

a. 

1V35  i*iy 

'“-.54  912*33' 

"o'. 

2 

1 

w  • 

_ .15312  *C3 

_ ; i  34  91 S  * 03_ 

_ - . 2C56E*Gi 

3 

1 

C< 

.  15672*03' 

-•  5491  a  *'<i  3 

-• 4323a  *Ct 

4 

_1 

_ 6  # _ _  _ _ 

_  •  16 33  F  *C  3 

_  5491  £  *33 

- •  5 925E  *C  1 

5 

i 

C  • 

.  16  99a  *C3 

“• 549lcf.3 

-•  7  7 56E  *ii  1 

_ 6 _ 

a. _ 

•  1/S5E  *C3 

-.3491c*:3 

95292*0  1 

7 

i 

a  • 

.  HS9E*G3 

-•  5'*9lt*-3 

- • 1221a  *0  2 

S 

_L 

0. 

•  19722  *C3 

54912*33 

1.473  F  *u  2 
-•  1725c  *v»2 

9 

l 

G. 

•  2a  756  *C3 

-.549if*i3 

.10 _ 

.  l 

a  • _ 

•  21792*03 

-.54912*03 

-.19652*02 

11 

l 

o  • 

.22R22*C3 

-•  5*»9l  E  *33 

-•  2198b  *02 

12 

i 

.  o* 

•  2366c  *C  3 

-•549lE*-:3 

-. 2426E*C2 

13 

1 

0  • 

.24392*03 

-.549l£*33 

-.26432*02 

14 

l 

a. 

•  25 97a  *G 3 

_-_._54FlE_*C3_ 

—  • 2  365c *02 

IS 

l 

a. 

•  2696  a  *C  3~ 

*• 54912*03 

-TOO  752  *0  2 

16 

i 

0. 

•23CCP*C3 

5491E*vJ3 

-.32382*02 

1 

2 

•3963&+0? 

•  14T5t*G3 

-• s357c  *33 

0. 

2 

2 

.3663*’ *02 

.15JIEK3 

-•  5357E*C3 

-. 2u35E*G1 

3 

2 

•  336U*C2 

•  15  o7c  *G3 

-•  53  5VE+G3 

-.47722*01 

SR  (I,J» 


•  25  7  8E  *  0  0 

.  js:2e*:q 


.34912.- 

.3**5G2* 

.3*172* 

.J37a6* 


_..T363P>0  2  .  1633  c  *0  3 -.  5T57E*C3 

•  3863c  *Q  2  .1699F*03  -•5357&*-03~ 


_;•  70 19c  *0 1 
-.91895*0  1~ 


•  3 3  3 It* 

•  33 19E*  0  0 

•  33 1  la  *  0  G 

•  33  3  82  *0  C 

•  33  0  7E*0  a" 
_^331_CE*_0  0 

•  33  l  uc  *  3Q~ 

•  33  2 1£*  0  0 

•  *3  5  8E  *  J  Q 

•  42  5  72*  G  0 
.4195E*00 

_.wl_J6E*0_G 
.  *0  9  8  F  *‘0  C” 


FIGURE  B-7.  PARTIAL  SAMPLE  OUTPUT  FROM  RAD FUR  FOR  THE 
INPUT  DATA  GIVEN  IN  FIGURES  B-5  and  3-6 
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O  O  Cj|c.  csl 


Tha  fall owing  thraa  data  card*  should  b«  includad  la  tha  Input  data 
strsaa  only  whoa  tho  first  varlabla  (A)  la  tha  laput  list  to  RFQRCE  Is  sat  to 
T&CE  (T). 


S.  Fifth  Data  Card.  Format  (3F10.4) 


Coluaas  1-10:  AA,  Dlaaatar  of  tha  cooling  plpa  la  aaadral,  bib  (30.) 
11-20:  AS,  Effaetlva  langth  of  tha  aaadral,  aa  (200.) 

21-30:  AC,  Laagth  of  tha  tubular  prafora,  aa  (200.) 

31-40:  AO,  Nuabar  of  forging  strokaa,  par  atfauta  (200.) 

41-30:  AE,  Fraetioa  of  stroka  whlla  tools  coatact  aatarlal  (0.S) 
51-40:  AF,  Ialcial  taaparatura  of  tha  dla,  C  (50.) 

61-70:  AC,  Zaltlal  taaparatura  of  tha  aaadral,  C  (50.) 

71-80:  AH,  Taaparatura  of  tha  aaadral  coolaat,  C  (50.) 


6.  Sixth  Oata  Card,  Forest  (8F1Q.4) 

Coluaas  1-10:  AK,  Aablaat  taaparatura,  C  (20.) 

11-20:  AL,  Seals  thiekaass  batvaaa  tha  blllat  and  aaadral,  aa  (0. 


7.  Savanth  Oata  Card.  Format  (8110)  -  -  - 

Coluaas  1-10:  / AfC"~Liaiting  nuabar  of  iterations  par  grid  movamanc  of 

(f  \gha  blllat  (20)  - ■*' 

11-20 1  ANx7^raB.« v r-i  H  '!i«<  vhara  siaulatloa  stops  ^9)  . 
21-30:  AO,  Nuabar  of  daslrad  Itaratloa  for  dabugglag  purposa. 


Laava  blank  otharvlsa. (0) 


31-40:  A?,  Nuabar  of  transvarsa  grid  Haas  la  tha  blllat  (12) 

41-50:  AQ,  Nuabar  of  longitudinal  grid  liaas  In  tha  oaadral  (5) 
31-60:  AH,  Nuabar  of  longitudinal  grid  liaas  la  tha  dlas  (4) 


Aa  aa  axaapla,  tha  input  data,  as  priatad  by  prograa  RFTEMF,  ara  shown 
In  Figura  B-8.  Tha  saapla  output  for  tha  saaa  input  data  is  givan  in  Flgura  3-9. 
Flgura  B-9a  is  a  partial  output  aftar  axtaaaiva  haat  traasfar  analysis.  (Whan 
tha  sacond  varlabla  B,  of  tha  first  card  is-  sat  "TRUE”) .  Figura  B-9b  is  tha 
output  obtainad  aftar  slsplifiad  haat  traasfar  aaalysis.  (Whan  tha  sacond 
varlabla  B  of  tha  first  card  is  sat  to  "FALSE"). 
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F ICDRIi  B-9a.  PARTIAI.  SAMPLE  OUTPUT  PROM  RFTEMP  POR  THE  INPUTS 
SHOWN  IN  PICHMES  R-5  AND  B-B  (Oetalled  Heat 
Transfer  Analysis) 


FIGURE  9b.  SAMPLE  OUTPUT  AFTER  SIMPLIFIED  MEAT  TRANSFER  ANALYSIS 


Th*  following  data  card*  chould  not  b«  included  in  eh*  input  data 
scream,  if  th*  -third  variable  (C)  in  the  input  list  to  RFORGE  is  set  to  FALSE  (F). 

8.  Eighth  Data  Card.  Foraat  (IIP) 

Colusa*  1-10:  AS,  Humber  of  inlet  angles  on  the  die  face  (3) 

9.  Ninth  Data  Card,  Foraat  (8F1Q.4) 

Column*  1-80:  AT,  Die  entrance  angle*.  Each  angle,  beginning  from 

th*  entrance  aide,  within  ten  columns,  degrees 
(20..  15.,  5.) 

Note:  Exclude  the  die  land.  Continue  on  an 
additional  card,  if  AS  >  8. 


10.  Tenth  Data  Card,  Format  (8F10.4) 

Column*  1-90:  AU,  Length  of  each  angular  segment.  Each  length 

beginning  from  th*  entrance  side,  within  ten 
columns,  ms  (150.,  225.,  30.) 

Note:  Exclude  the  die  land.  Continue  on  an 
additional  card,  if  AS  >  8. 

Aa  an  example,  the  input  data,  as  printed  by  the  program  RJC0MP,  are 
given  in  Figure  B-10.  Th*  sample  output  for  th*  same  input  data  is  given  in 
Figure  B-ll. 


DESCRIPTION  OF  RADFLW  ANT)  ASSOCIATED  SUBPROGRAMS 


This  sec  of  computer  programs  determines  a  kinematically- admissible 
upper-bound  velocity  field,  strain  races,  strains,  and  an  upper  bound  on  the 
total  forging  load  in  radial  forging  of  gun  barrels  over  a  stationary  mandrel. 
In  these  programs,  th*  new  international  system  of  units,  "System*  Interna¬ 
tional"  (SI),  has  been  used.  A  general  model  of  the  radial  forging  of  tubes, 
as  shown  schematically  in  Figure  B-12,  is  considered.  Th*  model  consists  of 
the  forging  and  of  the  sizing  zones.  Th*  upper-bound  method  was  utilized  for 
analyzing  the  deformation  process.  The  details  of  the  analysis  and  the  mathe¬ 
matical  derivations  are  given  in  Appendix  D  attached  to  the  final  report  of 
Phase- I  studies. 
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NUMBER  of  the  forging  TOOLS  *  4 

FORGING  TOOL  ENTRANCE  ANGLES .DEGREES 

20 .0 0  15.00  5.00  _ 
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nn.nn  in.Qfl _ _ _ _ _ 


FIGURE  B-10. 


SAMPLE  INPUT  DATA  AS  PRINTED  BY  RFCOMP 
PRIOR  TO  PERFORMING  COMPUTATIONS 
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FIGURE  B-I2.  SCHEMATIC  DIAGRAM  OF  RADIAL  FORGING  OF 
GUM  BARBELS  (WITHOUT  SINKING  ZONE) 


-Input  co  and  Oucout  from  Program  RADFLV 

The  computer  program  RADFLV  uses  as  input  chs  orocsss  variables 
raad  by  RFORGE.  In  addition,  RADFLV  raads  in  the  data  required  for  grid 
formation. 

The  computer  prograa  RADFLV  coapuces  the  following  quantities: 

(1)  Length  of  the  forging  zone 

(2)  Distance  of  the  neutral  plane  from  the  die  entrance 

(3)  Average  strain  and  average  strain  race 

(4)  Average  flow  stress 

(5)  Total  radial  forging  load  per  cool 

(6)  Energy  rate  supplied  by  a  forging  tool 

(7)  The  velocity  and  strain-race  distribution. 

The  prograa  RADFLV  is  also  applicable  co  the  radial  forging  of  rounds.  For  this 
purpose,  it  is  necessary  to  siaply  specify  the  diameter  of  the  mandrel  and  the 
friction  at  the  cube  aandrel  interface  as  zero. 

Details  of  RADFLV  and  Ocher  Subprograms 

The  basic  functions  of  the  rouclne  RADFLV  and  the  associated  subpro¬ 
grams  are  briefly  described  below. 

Subroutine  RADFLV 

This  subroutine  serves  essentially  as  a  coordinating  routine.  It 
collects  the  input  data,  calls  ocher  subprograas,  performs  some  simple  arith¬ 
metic  conversions,  sad  prints  out  the  final  resules  of  the  upper-bound  analysis. 


Function  ENERGY 

Purpose :  For  a  given  value  of  the  distance  (zq)  of  che 
neutral  plane  from  che  die  entrance,  this  sub¬ 
program  determines  che  total  energy  race  of  the 
process.  This  consists  of  che  power  of  plastic 
deformation,  che  power  co  overcome  frictional 
constraints,  and  che  shear  work  race  at  the 

boundaries  of  veloeicv  discontinuities. 
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Function  Reference:  POWER  •  ENERGY  (ZN) 

Input:  ZN  Is  ch«  distinct  of  the  neutral  plan*  from  -he  die 
entrance. 

Variables  Transferred  through  COMMON*:  ALFA,  AL1,  AL2,  CFRICD,  CFRICM,  DVEL , 

PAI,  RE,  RM,  RN ,  RO ,  TALFA 
Called  3v:  SIMPLE 

Sub rou tint  GSYST 

Purpose :  This  subroutine  establishes  an  orthogonal  rectangular 

grid  system  over  an  axial  cross  section.  GR(I)  and  GZCJ) 
give  the  radial  and  the  axial  grid  lines,  and  VOL(I,J)  are 
the  volumes  of  each  ring  element.  The  elements  on  the 
tapered  portion  of  the  tool  are  of  -triangular  cross  section. 

Calling  Sequence:  CALL  GSYST 

Variables  Transferred  through  COMMON:  GR,  GZ,  I END,  JEND,  NR,  NREXT,  NZ,  PAI, 

RE,  RM,  R2f,  RO,  VOL 
RADFLW 


This  subroutine  determines  the  optimum  value  of  the  unknown 
parameter  zq  (distance  of  the  neutral  plane  from  entrance) 
by  minimizing  the  total  energy  race  of  the  process.  A  sim¬ 
plex  method  for  function  minimization  is  used. 

CALL  SIMPLE  (X,N,AL?A,3ETA, GAMA, XX, ELIMIT, LIMIT, NPR) 

H  •  the  number  of  independent  variables 

X  ■  (N+l)  secs  of  initial  guesses  of  independent  variables 
ALFA  •  reflection  coefficient,  a  positive  constant,  less  chan  1 
BETA  *  contraction  coefficient,  a  positive  constant,  less  chan  1 
GAMA  ■  expansion  coefficient,  a  positive  constant,  greater  chan  1 
EXIMIT  -  error  limit  on  output  XX 
LIMIT  •  a  limit  on  the  number  of  iterations 

NPR  ■  a  printing  code  during  minimization,  1  for  yes,  0  for  no. 

*  Variables  la  the  COMMON  are  described  at  the  end  of  this  appendix. 


Called  by: 

Subroutine  SIMPLE 

Purpose: 

Calling  Sequence: 

Input: 
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Output:  XX  •  an  Array  of  optimum  independent  variables. 

Cal lint  RoueinAa;  ENERGY 
Callad  Bv:  RAC FLU 

Subrouting  MAI ERL 

Purposa :  This  aubroueina  furnishes  cha  eabulacad  flov  strass 

data  for  a  ranga  of  strains,  strain  ratas,  and  tetnpera- 
turas.  It  also  ealculatas  tha  flov  strass  for  given 
values  of  strain,  strain  rata,  and  temperature. 

Calling  Saouanea ;  CALL  MATERL  (STRAIN,  STRSAT,  TEMP,  FSTRES) 

Input :  STRAIN  •  strain,  SIR RAT  •  strain  rata,  TEMP  •  taoparatura 
Output :  FSTRES  •  flow  strass 

Variables  Transferred  through  COMMON:  DENSTY,  SPHEAT 
Calling  Routines:  FSTRS 

Callad  Bv:  FLOW,  RaDFLW,  RFCOMP,  RFOR LD,  TDEFOM 

Function  FSTRS 

Purposa:  Using  the  flow-stress  data  available  from  the  subroutine 
MAIERL,  this  subprograa  performs  cha  interpolation  for 
determining  the  flow  stress  for  given  values  of  strain, 
strain  race,  and  temperature,  which  arc  also  supplied  by 
the  subroutine  MAIERL. 

Function  Reference:  FSTRES  ■  FSTRS  (MSTRN,  MTEMP,  T,  ASTR,  AC,  AM,  TEMP,  STR, 

STRSAT) 

Input:  TEMP  •  temperature,  STR  ■  strain,  STRSAT  »  strain  race. 

The  ocher  variables  MSTRN,  MTEMP,  T,  ASTR,  AC,  and  AM  are 
internally  defined  within  cha  subroutine  MATERL. 

Calling  Routines :  AITO 
Callad  By:  MAIERL 

Subroutine  AITXN 


Purpose:  This  subroutine  performs  polynomial  interpolation  of  a 

cable  of  values  7  versus  X. 
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CALL  AITXN  (X.  Y ,  N,  X,  XB ,  YB ,  I,  IEX) 

X  ■  a  one-dimensional  array  of  monotonous  independent 
variab la 

Y  ■  a  one-dimensional  array  of  dependent  variables 
X  ■  cha  auabar  of  X-Y  pairs 

X  •  cha  dagraa  of  Lacarpolaclag  polynomial,  X  <  X  -  1 
XB  •  cha  valua  of  cha  ladapaadaac  variable  co  be  used 
Oucpuc :  YB  ■  cha  lacarpolacad  valua  of  cha  dapandaac  variable 
IEX  ■  iacagar  for  error  cheek,  1  if  extrapolation  was 
performed,  0  otherwise. 

The  variable  naaa  T  raprasaacs  a  one-dimensional  array 
of  ac  lease  2  (K+l)  words  used  for  camporary  storage. 
Called  Bv:  FSTRS ,  RFCOMP 


Calling  Sequence : 

Input; 


DESCRIPTION  OF  RrTSy?  AND  ASSOCIATES  SUBPROGRAMS 


RFTEXP  and  cha  associacad  rouciaas  are  used  co  obcain  temper acure 
discribu cions  la  cha  billec,  cha  mandrel,  cha  die,  and  cha  produce  la  radial 
forging  of  tubes.  The  temperature  analysis  used  in  developing  these  programs 
utilizes  a  fiaica-diff arenca  algorithm.  The  numerical  procedure  approximates 
the  heat  generation  and  the  heat  transportation  during  a  time  interval  taking 
place  instantaneously  followed  by  cha  simultaneous  heat  conduction  during  cha 
same  time  interval.  The  repetition  of  chase  two  steps  simulates  numerically 
cha  deformation  process  and  gives  the  temperature  distribution  as  a  function 
of  time .  The  detailed  derivations  of  the  mathematical  equations  used  in  these 
programs  are  given  in  Appendix  B,  attached  co  the  final  report  of  Phase-II 
studies.  A  brief  description  for  RTTEXP  and  its  associacad  subroutines  are 
included  later  in  this  appendix. 
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Input  to  and  Output  from  RJTSyP 


The  process  variables  read  by  RFORGE  arc  usad  by  RFTEM.  All  the 
velocity-related  data  gaaaratad  by  RADFLW  art  raad  froa  a  storage  flit.  In 
addition,  RFTEyp  raada  la: 

(a)  Initial  unifora  taaparaturaa  la  the  dla,  tha  aandral, 
tha  coolant  and  tha  ataosphara 

(b)  Billat  langth  and  aandral  langth 

(c)  Number  of  grid  linas  in  tha  billat,  dia,  and  aandral 

(d)  Thlcknaa*  of  tha  scala  baevaan  aandral  and  tuba 

(a)  Fraction  of  tha  stroka  during  which  tool  kaaps  contact 
with  tha  aatarial. 

Tha  printout  froa  RFTEMP  is  in  tha  fora  of  tablas  which  lncluda 
variables  as  raad  in,  tha  velocity-related  data  and  temperatures  in  the  billat, 
aandral  and  dia  at  various  else  intervals.  Tha  temperatures  are  also  printed 
at  the  and  of  tha  simulation. 


Details  of  RFTEMP  and  Other  Subprogram* 


Tha  basic  functions  of  RFTDfP  and  tha  associated  subprograms  are 
described  below. 

Subroutine  RFTEMP 


In  addition  to  monitoring  the  execution  of  its  associated  routines , 
this  subroutine  reads  tha  process  variables,  prints  all  tha  input  information 
with  appropriate  captions  for  tha  purpose  of  verification,  and  reads  tha 
velocity-related  data.  It  modifies  the  velocity  field  to  suit  tha  units  in 
temperature  prograa  and  prints  with  suitable  captions. 

Subroutine  COEFF 


Purpose:  This  subroutine  calculates  the  necessary  parameter 
values  for  heat-transfer  calculations. 


Calling  Sequence:  CALL  CO  EFT  (So  formal  parameters) 

Variables  Transferred  through  COMMON:  AMBTEM,  BILLES,  BILTSM,  CLZAR2 ,  COLTEM, 

OELIAZ,  DELZ,  OIETEM,  FRAC,  GRX,  I3DIE,  EEDIE,  MANLO, 

MAST EM,  MAXL,  MAXLM1 ,  MAXLP1,  MDIE,  MPRD,  MPRDM1,  MPRDP1, 
MTOT,  MTOTM1,  MTOT? I ,  SBC.,  MDIE,  MDISM1,  NDIEP1,  MMAN, 
MMANM1,  MMANP1,  MMANT,  MRDL,  MTOT,  MTOTP1,  HZ,  MZIEMP. 
RTEMP,  SI,  R2,  R21,  R3 ,  R32,  R4 ,  R43,  SLUT,  T,  TSCR  and  Z. 
Cal lad  By:  RFTEMP 

Subroutlna  FORGSG 


Purpose:  Based  on  stability  criterion,  chia  subroutine  estimates 

cha  critical  tiaa  interval,  it,  for  compiling  temperatures 
by  tha  finite-difference  method.  In  addition,  FORGSG  cal¬ 
culates  such  quantities  as  the  Instantaneous  billet  volume, 
volume  of  mecal  trapped  in  tha  deformation  zone,  product 
volume,  and  the  number  of  steps  required  to  move  the  billet 
by  one  grid  length. 


Calling  Sequence:  CALL  FORGNG 

Variables  Transferred  through  COMMON;  AXA ,  AKB,  AKD,  Aid,  3ILIR,  BILOR,  BILVEL, 

CANG,  CLEAR2 ,  COLTEM,  CONST,  CONVOL,  CHA,  CRB,  CRD,  CRM, 
DELT,  DELTAZ,  DELZ,  DIETEM,  FRAC,  GRX,  GRY,  HTCB ,  HTCD , 
HTFD,  IB) IE,  IPRI3T,  ITESIM,  JTEMP,  KJ,  LIMIT,  MAXL, 

MBPROD,  MDIP,  MEPROD,  MNPRDO,  MNPROD,  MPRD,  MTOT,  S3BIL,  MB 
NCONT,  MDIE,  MDIDU,  MDIEP1 ,  MDOASM,  MEBIL,  MMAN,  MMANP1, 
NKANT,  MRDL,  MSTRK,  STOT,  MZTEMP,  IP,  POWER,  RTEMP,  Rl,  R21 
R32 ,  R4 ,  R43,  SURT,  T,  X7RAC,  TSCR,  TCBIR,  TUBOR,  VR,  VZ,  Z 
Calling  Routines:  AIR,  BILLET,  DIE,  FRGTEM,  MANDL,  RFORLD,  TDEFQM,  TPRINT 
Called  By:  RFTEMP 

Subroutine  FRGTEM 

Purpose:  FRGTEM  estimates  the  temperature  distribution  after  heat 
transfer  during  radial-forging  operation. 
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Calling  Sequence :  CALL  FRGTEM 

Variables  Transferred  through  COMMON:  AKA ,  AKB ,  AKD,  A KM,  AMBTEM,  CLEAR2 ,  CLEARS  , 


COLTEM,  CRA,  CM,  CRD,  CRM,  DELT,  DELTAZ,  DELZ,  FRAC ,  HTCB , 
HTCD,  HTCM,  HTFB ,  HTFD ,  IBDIE,  ISDIE,  IPRINT,  KJ,  MAXL, 
MEPROD,  MNP  ROD ,  MPRD ,  MTOT,  NBBIL,  NCONT,  NDIE,  NEBIL, 

NMANT,  NTOT ,  NZTEMP ,  RTSMP ,  R1 ,  R21 ,  R32,  RA ,  R42,  R43. 

SURT,  T,  TSCR 

Calling  Routings: 

AIR,  BILLET.  DIE,  MANDL 

Called  By: 

FORGNG 

Subroutine  SIMPHT 

Purpose : 

SIMPHT  estimates  by  simplified  approximate  equations,  (a) 

average  temperature  of  the  deformation  zone  at  the  end  of 

the  process,  (b)  heat  gained  due  to  metal  deformation,  (c) 

heat  lost  to  atmosphere  by  convection,  (d)  heat  lost  to 

atmosphere  by  radiation,  and  (e)  heat  lost  to  mandrel  by 

heat  conduction. 

Calling  Secuence: 

CALL  SIMPHT 

Variables  Transferred  through  COMMON:  AKD.  AMBTEM.  BILOR.  BELTEM.  BILVEL .  CCLTEM 

CM,  01,  DLAND,  DO,  MANTEL,  PI,  POWER,  PRESSR,  TANG,  TUBIR, 

TUBOR 

Calling  Routines: 

BILLET,  MANDL 

Called  By: 

RITEMP 

Subroutine  TDEFOM 

Purpose : 

TDEFOM  calculates  the  temperature  gain  due  to  plastic 

deformation  of  the  billet  material. 

Calling  Sequence: 

CALL  TDEFOM 

Variables  Transferred  through  COMMON:  AMBTEM.  BILOR.  BILVEL.  CFRICD.  CFRICM. 

CM,  CRD,  CRM,  DELT,  DELIAZ,  DELZ,  DLAND,  ENVLOP,  GRX, 

CRT,  IBDIE,  I ED IE,  IPRINT,  KJ,  MAXL,  MEPROD,  MNPROD ,  MPRD, 

MTOT,  NBBH ,  NCONT,  NDIE,  NEBIL,  NMANT,  NRDL,  NZIEMP,  RTEHP, 

STRAIN,  STRRAl,  SURT,  T,  TSCR,  TUBIR,  TUBOR,  TUBVEL,  VR,  and 

VZ  as 
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Calling  Routines:  BILLET,  0t£,  MAMDL,  MATERL 


Called  Bv:  FORGNG 

Subroutine  RFORLD 

Purpose :  Calculates  che  radial-forging  load  par  cool. 

Calling  Sequence:  CALL  RFORLD 

Variables  Transferred  through  COS*(ON:  ALOAD,  CFRICD,  CFRICM,  CLEAR2 ,  COMVOL, 

OELTAZ ,  DVEL,  MBBIL,  SID II,  NMANT,  NRDL,  MZTSMP,  PI,  POWER, 
RIEMP,  STRAIN ,  S TREAT,  TSCR,  TUB II,  TUBOR,  and  VZ 
Calling  Routines:  MATSRL 
Callad  By:  FORGNG 

Subroutine  AIR 

Purpose:  Defines  cha  thermal  properties  for  eha  atmospheric  air. 

Calling  Saquanca:  CALL  AIR (TEMP) 

TEff:  Temperature  of  cha  air 

Variables  Tranafarrad  chrough  COMMON:  AKA,  CRA 
Callad  By:  FORGNG,  FRGTEM 

Subroutine  BILLET 

Purpose :  Dafiaaa  eha  thermal  propartias  for  cha  billet  material. 
Calling  Sequence:  CALL  BILLET(TEMP) 

TEMP:  Average  camper acura  of  cha  billet 
Variables  Transferred  through  COMMON:  AKB ,  CRB ,  HTCS ,  3TF3 
Called  By:  FORGNG,  FRGTEM,  SD2PHT,  TDErOM 


Subroutine  DIE 


Purpose:  Defines  the  thermal  properties  for  the  die  material. 
Calling  Sequence:  CALL  DIE (TEMP) 

TEMP:  Average  temperature  of  the  die 
Variables  Transferred  through  COMMON:  AKD,  CRD,  HTCD,  HT7D 

Callad  By:  FORGNG,  FRGTEM,  TDEFOM 
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Subroutine  MANDL 


Purpose:  Defines  the  thermal  properties  for  the  mandrel  material. 
Calling  Sequence:  CALL  MANDL  (TEMP) 

TEMP:  Average  temperature  of  the  mandrel 
Variables  Transferred  through  COMMON :  Aid,  CRM,  HTCM,  HTFM 
Called  Bv:  FORGNG,  FORTEM,  SIMPHT,  TDEFOM 

In  addition,  RFTEMP  utilizes  subprograms  MATERL,  FSTRS ,  and  AITKN 
described  under  RADFLW. 


DESCRIPTION  OF  RFCOMP  AND  ASSOCIATES  SUBPROGRAMS 


The  computer  program  RFCOMP  determines  the  quantitative  relationships 
between  various  process  variables  in  radial  forging  of  gun  barrels,  or  tubes, 
using  the  compound-angle  dies  (multiple  angles  on  the  die  entrance) .  In  its 
present  structure,  RFCOMP  can  consider  up  to  nine  different  angles  on  the  die 
entrance  and  is  applicable  to  both  the  round-faced  and  the  flat-faced  forging 
cools.  A  general  model  of  radial  forging  of  tubes,  as  shown  in  Figure  B-13  with 
three  entrance  angles,  is  considered.  This  model  consists  of  the  forging  and  the 
sizing  zones.  The  slab  method  was  utilized  for  analyzing  the  deformation  meehani 
The  program  RFCOMP  is  also  applicable  to  the  radial  forging  of  rounds. 
For  this  purpose,  it  is  necessary  to  simply  specify  the  Inside  diameter  of  the 
billet,  the  diameter  of  the  mandrel,  and  the  friction  at  the  cube-mandrel  inter¬ 
face  as  zero. 


Input  to  and  Output  from  RFCOMP 

In  addition  to  the  process  variables  read  by  RFORGE,  RFCOMP  uses  as 

input: 

(a)  Number  of  angular  segments  in  the  dies 

(b)  Entrance  angles,  degree 

(c)  Lengths  of  angular  segments  in  dies,  mm 
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Sinking 


Tha  computer  program  RTCOMP  compueas  cha  following  quancicias: 

(a)  Langehs  of  cha  forging  and  eha  sizing  zonas 

(b)  Discanca  of  eha  naueral  plana  from  eha  dia  ancrance 

(c)  Serain,  serain  raca,  flow  serass,  axial  serass,  and 
eha  radial  prassura  diseribueion  along  eha  axis  in 
eha  daformaeion  zona 

(d)  Tocal  radial-forging  load 

(a)  Avaraga  radial  forging  prassura 

(f)  Componanes  of  eha  eoeal  load  in  eha  fora  of  forging 
and  sizing  loads 


Daealls  of  RTCOMP  and  eha  Assoeiaead  Subprogram 

Tha  basic  funccions  of  RFCO**?  and  cha  various  subprograms  of  RTCOMP 
ara  briafly  dascribad  balow. 

Sub roue ine  RFCOM? 


This  roucina  sarvas  assancially  as  a  coordinacing  roucine.  Ic 
collaces  cha  inpue  data,  calls  oeher  subprograms,  parforms  some  simple  arieh- 
macic  conversions,  and  princs  oue  eha  final  rasules  of  che  serass  and  load 
analysis . 


Sub roucina  PRESUR 


Purposa: 


Calling  Saouanca; 

Inpue: 


This  subroueina  ganaraeas  eha  radial  prassura  and  che 
axial  serass  diserlbueions  scarring  from  aichar  and  of 
eha  daformaeion  zona.  Tha  lnearsaceion  of  cha  evo  pras¬ 
sura  curvas ,  ehus  ganaraead,  glvas  eha  locaeion  of  cha 
naueral  plana  and  minimum  of  eha  evo  values  of  eha  cal- 
culaead  radial  prassura  is  cakan  as  eha  aeeual  prassura. 
CALL  PRESUR  (I CALL) 

I CALL  is  alehar  1  or  2.  Whan  ICALL  ■  1,  eha  subroueina 
PRESUR  is  eallad  in  ordar  eo  daearmina  eha  naueral  plane. 
Whan  ICALL  •  2,  ie  is  eallad  for  deearmining  eha  aeeual 

radial  prassura  diseribueion. 
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Variables  Transferred  chroujh  COWON:  ALFA,  BPULL,  CFRICD,  CFRICM,  3 LAND, 

FPULL,  FSTR,  FSTRES ,  NFORG,  NSTEP ,  PAL,  PAR,  R,  ri ,  R2 , 
STR,  STRAIN,  STAR,  STRRAT,  SZL ,  SZA,  TL,  Tl,  Z 
Calling  Routines:  FLOW 
Called  3v:  AFC CMP 


Subroucin*  PLOW 

Pumas* :  Ac  a  given  radial  plans  in  eh*  deformation  eon*,  this 
subroucin*  calculates  eh*  scrain,  eh*  strain  rat*,  and 
it  d*c*rsiin«s  eh*  flow  s cress  by  calling  eh*  subroucin* 

matesl: 

Calling  S«ou*nc*:  CALL  FLOW  (ZZ,  RR) 

Inout:  ZZ  is  eh*  distance  of  a  given  radial  plan*  from  eh*  di* 
entrance,  and  RR  is  th*  outside  radius  at  chat  location. 
Variables  Transferred  through  COMMON:  DLAND,  DVEL,  FSTRES,  RI,  RM,  RN,  RO, 

R2,  STRAIN,  STRRAT,  TSMPR,  TL 
Calling  Routin*:  MATESL 

Called  3v:  PRESUR,  RFCOMP  . 

In  addition,  RFCOMP  utilizes  the  subprograms  MAIZRL,  FSTRS ,  and 
AlTKN  described  under  RADFLW. 
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TABLI  8-1.  CROSS  REFERENCE  TABLE  FOR  SUBPROGRAM  CALLINGS 


Calling 

Calltd  By 

1.  AIR 

FORGNG , FRGTEM 

2 .  AlTKN 

FSTRS, RFCOMP 

3 .  BILLET 

FORGNG ,  FRGTE** ,  TDEFOM 

4 .  COEFF 

RFTEMP 

5.  DIE 

FORGNG , FRGTEM , TDEFOM 

6.  ENERGY 

SIMPLE 

7.  FLOW 

MATERL 

PRESUR, RFCOMP 

8.  FORGNG 

AIR. BILLET , DIE , FRCTEM , MANDL . 
RFORLD , TDEFOM, TPRINT 

RFTEMP 

9.  FRCTEM 

AIR .BILLET . DIE , MANDL 

FORGNG 

10.  FSTRS 

AlTKN 

MATERL 

11.  GSYST 

RADFLW 

12.  MANDL 

FORGNG , FRGTEM , TDEFOM 

13.  MAT ERL 

FSTRS 

FLOW ,  RADFLW , RFCOMP , RFORLD , 
TDEFOM 

14.  PRESUR 

FLOW 

RFCOMP 

15.  RADFLW 

GSYST ,  MATERL ,  S IMPLE 

RFORGE 

16.  RFCOMP 

AlTKN, FLOW, MATERL, PRESUR 

RFORGE 

17.  RFORGE 

RADFLW , RFCOMP , RFTEMP 

18.  RFORLD 

MATERL 

FORGNG 

19.  RFTEMP 

COEFF, FORGNG 

RFORGE 

20.  SIMPHT 

BILLET, MANDL 

RFTEMP 

21.  SIMPLE 

ENERGY 

RADFLW 

22.  TDEFOM 

BILLET, DIE, MANDL, MATERL 

FORGNG 

23.  TPRUTT 

FORGNG 
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TABLE  B-2.  LIST  OF  VARIABLES  PRINTED  FOR 
DIFFERENT  SETTINGS  OF  I PRINT 


The  following  variabl as  will  b«  printed  whan  IPRINT  is  sac  higher 
or  equal  co  cha  value  shown  balow. 


Talus  of  IPRINT 


Tarlablas  Listed 


0 


L 


(a)  All  arror  messages 

(b)  DO,DI,DE,DM,DVEL,AXFEED,DIlANG,3PULL, 
F7ULL , CFRI CD , CFRICM , 3 ILTEM , PHI , S  AI , 
and  ITOOL  from  RFORGE 

(c)  CF.DLAND.ZN.RN.ENVLOP,  and  PN  from 
RFCOMP 

(d)  2N.RN.S TRAIN, STRRAI , FSTRES , SN , PN , 

ALOAD .ANGPR , FORGL , PCFORG , S I2EL , PS I2E , 
and  POWER  from  RFCOMP 

(a)  Tamparatura  discrlbuclon  in  the  biiiec, 
Cha  dla,  che  mandrel  and  cha  produce  ae 
Cha  baginning,  ac  times  when  billac 
novas  chrough  one  grid  line  discance, 
and  ac  che  end  of  ieeraeions  from 
RFTEMP 

(a)  TELCTT ,TMPRTR , PRESSR , IPRINT  from  RFORGE 

(b)  NR.NZ.NRECT  from  RADFLW 

( c)  AL1 ,2N,  STRAIN  ,  STRRAT ,  FSTRES  ,  ALOAD ,  PRESL’R , 
TFIN,  and  TWORK  from  RADFLW 

(d  )  DP  IPE ,  MANLEN ,  3 ILLEN ,  STROKE ,  AMBTEM ,  D IETEM , 
MANTEM ,  COLTEM ,  NB ILT ,  NMANT ,  MD  IE ,  L IMIT , 
SCALE,  and  TFRAC  from  RFTEMP 

(a)  NANGLE  from  RFCOMP 
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2  (a)  ALFA(J) ,AL(J)  J  •  L,  N ANGLE  from  RFCO*? 

3  (a)  I , J,GZ( J) , GR(I) ,UZ(I , J) ,UR(I,J) ,SR(I,J) 

for  I  -  1  co  IENOCJ)  and  J  •  1  to  NZ?I 
from  RADFLW 

(b)  NRDL,MAXL,MZ(J)  J  -  1  co  >HOL  from 
RFTEMP 

(c)  J,Z(J),VZ(I.J),VR(I,J),S?TUi;T(I,J)  ,  and 
STRRAKI , J)  for  I  •  1  co  NZ(J)  and 

J  ■  1  to  NRDL  from  RFTEXP 

4  (a)  I.J.ICP, IDS, RRG.ICLF.CEQ, CONST,  and 

CHECK  for  I  »  1  co  MTOT  and  J  •  2  co 
NTOTMI  from  FORGNG 

(b)  K,I,T(I,1) ,T(I+1,1) ,T(I,MBPROD)  for 
K  ■  1  co  NDO  and  I  ■  1  co  NZ(K)  from 
FORGNG 

(c)  I,J,KC?.KLHtKRG,ia.F,CEQ,TZl,+Z2,rS,rj?, 
TDN.TRG.TLF.TSCRC.J)  ,T(I,J)  for  I  -  1 
Co  JfTOT  and  J  •  NBBIL  to  JEND  from 
FRGTEM 

(d)  TRISE  for  J  •  1  Co  NDXE  and  I  -  MMAN7 
co  NZ( J)  from  TOEFOM 
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Oescriaeion  of  the  Variables  in  the  CC-yyCN  Statements 


Most  of  the  important  variables  are  transferred  beeveen  subroutines 
through  COMMON  statements,  named  ARRY  and  SCLR.  For  easy  reference,  all  : r.e 
variables  in  the  COMMON  statements  are  described  below  in  alphabetical  crier . 


AXA 

AK3 

AKD 

AJOf 

AL 


All 

AL2 

ALFA 

ALFA 

ALOAD 
AMBTEM 
A R 


AXFEED 

3ILL53 

3IL0R 

BILTEM 

BIL7EL 

BPTJLL 

CANG 
CTRL  CD 
CFRICM 

a 

CLEAR2 


Thermal  conductivity  of  air 

Thermal  conductivity  of  billet  material 

Thermal  conductivity  of  die  material 

Thermal  conductivity  of  mandrel  material 

Lengths  of  angular  segments  on  the  die 

face  with  multiple  angles 

Length  L^  in  Figure  3-5 

Length  in  Figure  B-5 

Die  angle.  Angle  a  in  Figure  B-5 

Array  of  angles  (in  RFCCMP  and  associated 

routines) 

Forging  load  per  tool 
Ambient  temperature 

Outside  radius  at  the  beginning  of  each 
angular  segment 

Axial  feed  of  the  billet  per  stroke 
Billet  length 
Billet  outer  radius 
Billet  temperature 
Billet  velocity 

Back-pull  force  (with  a  negative  sign  if  back 
push) 

Cosine  of  the  die  angle 
Friction  shear  factor  at  the  die  surface 
Friction  shear  factor  at  the  mandrel  surface 
Axial  distances  of  starting  points  of  each 
angular  segment 

One-half  of  the  clearance  beeveen  billet 
inner  radius  and  tube  inner  radius 
One-fourth  of  the  clearance  beeveen  billet 

inner  radius  and  tube  inner  radius 
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CLEAR* 


COLIZM 

CONVOL 

cra 

CRB 

CRD 

CRM 

DE 

DELT 

DELTA2 

DEL2 


DEN  STY 
DI 

DIEANG 

DIETEM 

DLAND 

DM 

DO 

DVEL 

DZ 

ENVLOP 

FPTJLL 

FR 

me 


Temperature  of  che  coolant  in  the  mandrel 
Volumes  of  ring  elements  in  the  deformation 
zone 

Product  of  specific  weight  and  heat  capacity 
of  air 

Product  of  specific  weight  and  heat  capacity 
of  billet  material 

Produce  of  specific  weight  and  heat  capacity 
of  die  material 

Product  of  specific  weight  and  heat  capacity 
of  mandrel  material 

Outside  diameter  of  the  forged  product 
Delta  T.  Time  interval  it  for  heat  transfer 
calculations 

Delta  Z.  Axial  distance  between  two  successive 
radial  grid  lines 

Similar  to  DELTAZ  but  has  a  fixed  value  which 

is  the  axial  distance  between  wo  successive 

radial  grid  lines  in  the  billet  outside  the 

deformation  zone 

Density  of  the  material 

Inside  diameter  of  the  tubular  preform 

Die  angle 

Die  temperature 

Die  land  length 

Inside  diameter  of  the  forged  product  or 
diameter  of  che  mandrel 
Outside  diameter  of  the  tubular  preform 
Radial  velocity  of  hammers 
Delta  Z  in  each  angular  segment 
Envelopment  factor  of  the  forging  tool 
Front-pull  force  (with  a  negative  sign  if 
front  push) 

Strain  rates/area  at  neutral  plane 
Have  values  equal  to  R^/AR^ 
where  R^  -  [RTEMP(I)  +  HXEff(I-l)  1/2 
AR,  -  RTEMP(I)  -  RTEMP(I-l) 
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FSTR  Array  of  flow  s cresses 

FSTRES  Avaraga  flow  seress  In  eha  nacarial 

GR  Radial  coordinacaa  of  grid  poincs 

GRX  Axial  diacancaa  of  successive  radial  grid 

linaa  scarring  froa  cha  axie  of  the  dia 

CRY  GRY(J)  is  cha  surfaca  araa  ac  cha  die-b iliac 

incarfaca  ac  grid  Poinc  J 

HTC3  Haae  cransfar  coafficianc  ac  cha  cylindrical 

surfaca  of  cha  billae  nacarial 

ETCD  Haac  cransfar  coafficianc  ac  cha  cylindrical 

surfaca  of  cha  dia  naeerial 

HTCM  Haac  cransfar  coafficianc  ac  cha  cylindrical 

surfaca  of  cha  nandrel  nacarial 

HTF3  Haae  cransfar  coafficianc  ac  cha  flac  surfaca 

of  cha  billae  nacarial 

HT7D  Haac  cransfar  coafficianc  ac  cha  flac  surface 

of  cha  dia  nacarial 

HTTM  Haac  cransfar  coafficienc  ac  cha  flac  surface 

of  cha  nandral  nacarial 

IBDIE  Subscrlpc  of  cha  radial  grid  lina  vich  which 

dia  back  saccion  begins 

IEOXE  Subscrlpc  of  eha  radial  grid  line  with  which 

dia  back  saccion  ends.  Soch  IBDIE  and  IED1E 
rafar  eha  scoraga  lo carious. 

IDTD  Number  of  radial  divisions  ac  an  axial  locaclon 

ins Ida  eha  dias 

IPRINT  Incermediaca  resulcs  ara  baing  prlnced  whan 

IPRIHT  is  sac  a  valua  baewaan  0  co  5 .  ?efer 
Co  Table  B-2,  page  B-33. 

ITESIJf  Maxima  allowed  nunbar  of  learacions  during  necal 

daforaacion 

ITOOL  A  coda  nuaber  for  cype  of  eool  (1  for  round  and 

2  for  flac  face) 

JEND  Humber  of  axial  divisions  ac  a  radial  locaclon 

inside  eha  dies 

Has  a  valua  equal  co  IEDIE  +  1  and  is  used  by 
TPRIST  for  prineing  purposes 
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JTEMP 


KJ 

LIMIT 


MANLEN 

MANTEM 

MAXL 

MAXL  Ml 
MAXLP1 
MBPROD 

MDIE 

MDIFF 


An  index  co  Indicate  finish  shape  produce 
has  begun  coming  out 

Subscript  of  the  radial  grid  line  at  which 
whan  the  billet  end  reaches,  the  process 
simulation  is  eo  be  stopped.  For  example, 
if  MDIE  is  equal  to  10  and  MBILT  Is  equal 
to  20,  the  billet  end  will  be  at  (10  +  20  •) 
30th  radial  grid  line.  To  scop  the 
process  when  one-half  of  the  billet  length 
is  forged,  LIMIT  must  be  set  to  (10  ♦  10") 

20.  To  forge  the  complete  length  of  the 
billet  LIMIT  must  be  10.  To  forge  the  billet 
through  only  one  grid  length  distance  LIMIT 
must  be  equal  to  29 . 

Mandrel  length.  Note:  real  variable 
Mandrel  temperature.  Mote:  real  variable 
Number  of  axial  grid  lines  from  the  center  of 
the  mandrel  to  the  outer  surface  of  the  billet 
MAXL  minus  one 
MAXL  plus  one 

Subscript  of  the  radial  grid  line  at  which 
product  begins  to  emerge 

Number  of  axial  grid  lines  in  the  die  above 
the  billet  outer  surface 
Number  of  radial  grid  lines  by  which  the 
temperature  distribution  in  the  mandrel  is 
shifted  with  respect  to  the  temperacure 
distribution  in  the  billet  at  the  end  of  each 
cycle  in  order  co  compensate  for  the  relative 
velocity  between  the  billet  and  mandrel 
Subscript  of  the  radial  grid  line  with  which 
product  front  face  coincides.  (Similar  to  I3DI 
and  IEDIE,  MBPROD  and  MEPROD  refer  to  the 
storage  locations). 
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MEPROD 


MNPRDO 


WSOD 

MPRD 

MPRDM1 

MPRDP1 

MTOT 

XTOTX1 

JfTOTPl 

MANGLE 

MBBIL 


MBILP1 

MBILT 

MCONT 


NDIE 

NDIEHL 

NDIEP1 

HDIV 

MDOASH 


MEBIL 


MFORG 

NMAN 

MMAtna 

MMANPl 


MNPROD  old.  The  value  of  cha  MNPROD  In  che 
previous  cycle  Is  scored  in  chis  variable 
Difference  becveen  ME? ROD  and  MBPROD 
Number  of  axial  grid  linos  in  che  produc: 

MPRO  minus  on* 

MPRD  plus  on* 

local  number  of  axial  grid  linas  including 
cha  dia 

MTOT  minus  ona 
MTOT  plus  ona 

Number  of  inlac  anglas  on  cha  dia  faca 
Subscripc  of  cha  radial  grid  lina  vich  which 
cha  b iliac  fronc  faca  coincidas  ac  cha  beginn- 
ing  of  a  cycla 
NBILT  plus  ona 

Numbar  of  radial  grid  lines  in  che  b iliac 
Number  of  concaccs  bacwaan  cha  die  and  che 
macarial  during  cha  calculaced  else  incerval, 

Ac 

Subscripc  of  che  radial  grid  line  vich  vhich 
cha  encrance  faca  of  cha  dia  coincide 
NDIE  minus  ona 
NDIE  plus  ona 

Number  of  divisions  in  each  angular  segmenc  (RFCOMP) 

MDO  assigned.  MDO  is  cha  nuabar  of  cycles  used 
for  process  siaulaeion.  By  assigning  a  value 
Co  NDOASN  nuabar  of  cycles  used  for  process 
siaulaeion  can  be  reduced  co  a  prasac  value. 

Useful  for  prograa  debugging. 

Subscripc  of  che  radial  grid  line  vich  which  che 
billec  and  face  coincides  ac  che  beginning  of  a 
cycle 

Nuabar  of  divisions  in  che  forging  zone 
Number  of  radial  grid  lines  in  che  mandral 
NMANT  minus  one 
NMANT^glus  one 


/v£aJ 

NKANT 

NR 

NROL 


NREXT 

NSTEP 

NSTRK 

NTOT 

N2 


N2TEVP 

PAI 

PHI 


PI 

POWER 

PR 

PRESUR 

PRL 

PRR 

R 


RE 

REDRAT 

RELOAD 

PI 

RM 


KJu~U<  4 

Number  of  axial  grid  lints  in  che  mandrel 
Number  of  radial  divisions  on  che  entrance  side 
Number  of  radial  grid  lines  from  die  exit  cs 
billet  end  surface 

Number  of  radial  divisions  on  che  exit  side 
Total  number  of  divisions  in  che  2one  of 
deformation 

Number  of  strokes  in  che  calculated  time  interval., 
At 

Total  number  of  radial  grid  lines 

N2CJ)  is  the  number  of  axial  grid  lines  at  che 

die  surface  for  che  J  radial  grid  line.  N2  is 

read  in  and  then  modified  to  account  for  che  grid 

lines  in  che  mandrel 

Same  as  N2 

Same  as  PI 

Tool  contact  angle  (for  round  tools  only) 

Refer  to  Figure  B-2 
A  constant  equal  to  3.141592654 
Forging  power  per  cool 
Array  of  radial  pressures  on  che  dies 
Logical  variable  to  indicate  chat  pressure 
and  load  calculations  are  to  be  performed 
Array  of  pressures  from  the  die  entrance  side 
Array  of  pressures  from  the  die  exit  side 
Radial  distances  of  che  axial  grid  lines  as 
read.  The  same  array  is  used  later  to  print 
radial  distances  in  mm  units,  also  outside 
radius  at  a  given  axial  location. 

Outer  radius  of  the  forged  product 
Area  reduction  ratio 
Not  used 

Inner  radius  of  the  tubular  preform 

Inner  radius  of  che  forged  product  or  radius 

of  the  mandrel 
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SN 

ao 


RPIPE 


RTEMP 


R1 

R 2 


R21 

R3 

R32 

R4 

R42 

R43 

SAX 

SCR 

SPHEAI 

SQ3 

SR 

SIR 

STRAIN 

STROKE 

STRX 

STRRAT 

SURT 


SZ 

SZI 

SZR 


Oitside  radius  ac  naucral  plana 

Outar  radius  of  the  tubular  prefora 

Radius  of  cha  cooling  plpa 

Modified  radial  distances  of  the  axial 

grid  lines  which  accounts  fir  cha  aandrel 

radius 

A  fixed  value  equal  to  (RTEMP(MMAN?i)  +■ 

RTEMP  (  NMANM1 ) )  /  2 

R1  minus  CLEAR2  (in  RFTEMP  and  associated 
subprograms) .  Same  as  RE  (in  RfCOMP  and 
associated  subprograms). 

(R2)2  -  (Rl)2 
R1  plus  CLEAR2 
(R3) 2  -  (R2) 2 

( RTEMP (NMANP1)  *  RTEMP (NMANT) ) /2 
(R4) 2  -  (R2)2 
(R4)  2  -  (R3) 2 

Tool  cangency  angle  (for  flat  cools  onlv) 

Scratch  array  for  determining  AR 
Specific  heat  of  cha  material 
A  constant  equal  to  cha  square  root  of  3 
Strain  rates  at  grid  points 
Array  of  strains 

Strain  in  cha  material  at  each  grid  point 
Number  of  forging  strokes  per  minute 
Array  of  strain  rates 

Strain  rate  in  the  material  at  each  grid  point 
Surface  area  of  the  element  at  the  intersection 
of  Ith  axial  grid  line  and  Jth  radial  grid  line. 
It  is  equal  to  [RTEMP(I)  >  RTEMP  (I  -  1)].'[2  x 
(RTEMP (I)  -  RTEMP (I  -  1))J 
Array  of  axial  stresses 

Array  of  axial  seresses  from  the  entrance  side 
Array  of  axial  stresses  from  the  exit  side 
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T 

TALFA 

TALFA 

TAMG 

TEM 

TEMP 

TERM 

TFRAC 

TL 

TMPRTR 

TO 

TSCR 


TUBIR 

TUBOR 

TUBVEL 

UR 

uz 

VELCTY 

VOL 

VR 

VZ 

Z 


zn 


Temperature  ac  ch«  intersection  of  ehe  Ith 

&xial  and  Jth  radial  grid  Unas 

Tangent  of  entrance  angle 

Array  of  tangents  of  angles  (in  RFCOM?  and 

associated  routines) 

Tangent  of  the  die  angle 
Array  of  temperatures 
Average  temperature 

Value  of  energy  rate  integral  in  RADFLV 
Fraction  of  the  stroke  time  during  which  cool 
and  material  are  in  contact 
Total  die  length,  equal  to 
Figure  B-10 

Logical  variable  to  indicate  that  temperature 
distribution  is  to  be  calculated 
RO  -  RI 

A  scratch  array  used  to  store  intermediate 

temperature  values 

Tube  inner  radius 

Tube  outer  radius 

Tube  velocity 

Radial  velocities  at  grid  points 
Axial  velocities  at  grid  points 
Logical  variable  to  indicate  that  velocity 
field  is  to  be  generated 
Volume  of  elements  around  grid  points 
Radial  velocity  of  the  element  at  the  intersec¬ 
tion  of  Ith  axial  and  Jth  radial  grid  lines 
Axial  velocity  of  the  element  at  the  intersection 
of  Ith  axial  and  Jth  radial  grid  lines 
Axial  distances  of  the  radial  grid  lines  as  read. 
The  same  array  is  used  later  to  print  axial 
distances  in  ma  units. 

Distance  of  the  neutral  plane  from  die  entrance. 
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